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ABSTRACT 

The  background  of  microstructure  and  its  effect  on  sound  propagation 
is  reviewed.   The  expendable  bathythermograph  (XBT)  system  and  current 
digitizing  procedures  are  evaluated  to  determine  their  adequacy  in 
investigating  vertical  thermal  microstructure .   Several  methods  of 
microthermal  analysis  are  proposed.   Temperature  gradient  versus  depth 
plots  are  used  to  analyze  data  from  several  ocean  regions  -  Atlantic 
Coastal,  Andaman  Sea,  Eastern  North  Pacific,  and  the  northern  boundary 
of  the  Gulf  Stream.   Analysis  of  these  data  suggest  a  water  mass  depend- 
ence of  the  thermal  microstructure.   The  significance  of  the  small-scale 
structure  to  Naval  operations  is  discussed. 
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I.   SIGNIFICANCE  OF  MICROTHERMAL  STRUCTURE 

Naval  interest  in  the  thermal  structure  of  the  ocean  is  presently 
centered  on  gross  features  such  as  sea  surface  temperature,  mixed  layer 
depth,  thermocline  gradient,  and  depth  to  the  bottom  of  the  thermocline. 
Some  attention  has  been  given,  especially  in  the  Canadian  Navy,  to  the 
transients  in  the  mixed  layer  [Tully,  1964]  and  their  possible  associ- 
ation with  internal  waves.   The  above  gross  features  are  used  primarily 
as  a  means  of  classifying  thermal  structure. 

Technological  developments  have  led  to  measurements  which  indicate 
the  thermal  structure  is  made-up  of  layers  of  varying  thickness;  from 
less  than  a  centimeter  to  tens  of  meters  [Stommel  and  Federov,  1967; 
Osborn,  1969;  for  example]. 

Of  the  variables  affecting  sound  velocity  in  the  ocean,  temperature 
has  the  greatest  effect.   The  presence  of  variations  in  temperature,  or 
temperature  gradients,  causes  sound  velocity  fluctuations  to  occur. 
Research  on  the  effects  of  the  thermal  layers  on  sound  in  the  ocean  show 
that  further  investigation  in  theory  and  experiment  is  required.   Small 
scale  sound  velocity  variations  cause  irregular  bending  of  sound  rays. 
The  microstructure  can  affect  attenuation,  forward  and  backward  scatter- 
ing, reverberation,  reflection,  and  refraction  which  can  converge  or 
diverge  the  energy  in  a  sound  beam  [Meyer  and  Romberg,  1963].   The 
physical  features  of  the  microthermal  structure  must  be  adequately 
described  before  their  effects  upon  acoustic  propagation  can  be  fully 
understood . 


The  purpose  of  this  thesis  is  to  review  the  present  knowledge  con- 
cerning microstructure  as  measured  by  expendable  bathythermographs  (XBT) , 
to  evaluate  the  capability  of  the  XBT  in  measuring  microthermal  structure, 
to  develop  a  scheme  of  quantitatively  examining  the  microthermal  structure, 
and  to  test  this  scheme  on  real  data. 


II.   REVIEW  OF  PRESENT  KN0WIJ5DGE 

A.   MICROTHERMAL  STRUCTURE  IN  THE  OCEAN 

Prior  to  1937,  vertical  thermal  structure  was  measured  at  discrete 
points  using  reversing  thermometers.   This  method  provided  a  point-to- 
point  analysis  of  the  vertical  structure  at  one  location.   Even  using 
point-to-point  measurements,  it  was  clear  that  the  temperature  vari- 
ability decreased  with  depth  in  most  regions.   In  fact,  the  most  common 
analysis  of  the  temperature  and  salinity  data,  T-S  correlation,  usually 
ignored  the  upper  100  meters  or  so.   However,  this  same  analysis  showed 
that  in  many  ocean  regions,  in  deep  water,  (e.g.,  below  1000  to  2000  m) , 
the  variability  was  on  the  order  of  the  instrument  error  [Sverdrup, 
Johnson  and  Flemming,  1942].   This  result  was  not  unexpected  because 
the  gradients  in  the  deep  ocean  became  smaller.   The  discrete  measure- 
ments of  the  Nansen  cast  are  not  suitable  for  studying  small  scale 
variability . 

The  development  of  the  mechanical  bathythermograph  (MBT)  by  Spilhaus 
[1938]  provided  the  first  convenient  means  of  measuring  a  continuous 
profile  of  temperature.   The  mechanical  bathythermograph  was  widely  used 
during  World  War  II.   A  large  number  of  bathythermograms  with  double 
traces,  which  could  not  be  explained  as  instrument  errors,  were  dis- 
covered [Urick  and  Searfoss ,  1948].   The  double  traces  indicated  that 
small-scale  horizontal  and/or  vertical  temperature  variations  occurred 
between  the  lowering  and  raising  of  the  MBT.   Urick  and  Searfoss  [1949] 
recognized  that  this  small-scale  microstructure  could  be  significant  to 
sound  propagation,  depending  on  the  size  and  shape  of  the  inhomogeneity . 


Horizontal  micro thermal  structure  was  measured  by  Holder  in  1944 
using  a  thermopile  mounted  on  a  submarine.   He  detected  horizontal 
temperature  differences  of  0.02F  over  distances  of  about  10  m  [Smith, 
1967].  Measurements  of  these  horizontal  temperature  variations  were 
also  made  by  Urick  and  Searfoss  [1948,  1949]  in  July  of  1948  near  Key 
West,  Florida.   With  a  thermocouple  mounted  on  the  conning  tower  and, 
later,  its  bow,  a  submarine  proceeded  at  a  constant  speed  and  depth, 
while  the  probe  measured  the  water  temperature  to  a  resolution  of 
0.003C. 
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Figure  1.   Microthermal  Variations  at  Various  Depths 
[Urick  and  Searfoss,  1948]. 


Urick  and  Searfoss  [1949]  calculated  the  peak-to-peak  root-mean- 
squara  (rms)  variation  for  peaks  in  excess  of  20  percent  of  the  total 
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variation  in  a  record  calling  this  the  root-mean-square  maximum  (rmsm) . 
A  significant  difference  between  the  value  of  rmsm  in  the  mixed  layer 
and  in  the  thermocline  was  observed.   In  the  mixed  layer  a  mean  rmsm 
value  of  0.013C  was  observed.   In  the  thermocline,  the  mean  value  was 
0.450C.   The  rmsm  increased  with  depth  in  the  mixed  layer  and  decreased 
with  depth  in  the  thermocline.   Applying  Taylor's  hypothesis,  the  mean 
horizontal  thermal  "patch  size"  was  found  to  vary  from  5  to  30  m  in  the 
mixed  layer  and  from  30  to  200  m  in  the  thermocline  region  [Urick  and 
Searfoss,  1949]. 

Liebermann  [1951]  mounted  a  sensitive  platinum-resistance  thermom- 
eter and  a  thermocouple  on  the  periscope  of  a  submarine.   The  platinum- 
resistance  element  could  resolve  spatial  temperature  variations  on  the 
order  of  10  cm  in  horizontal  distances  as  the  submarine  proceeded  at  a 
constant  speed  of  a  few  meters  per  second  at  a  fixed  depth  (Figure  2). 
To  characterize  the  temperature  fluctuations,  he  used  the  autocorrelation 
function.   "Patch  size"  was  defined  as  the  distance  of  the  microthermal 
structure,  within  which  the  temperature  variations  maintained  coherence. 

The  function  is  given  by: 

CD 

R  (d)  =  J   [T  (x)-T  (x+pQ)]  dx.  (1) 

0 

His    results    could   be    fit   equally  well  by   an  exponential, 

R    (P)    =  e"p/a  (2) 

or   the   Gaussian, 

.2 


r    (D)   =  e"(p/a)'  (3) 


where 
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T  (x)  =  temperatute  at  some  point  x; 
T  (x+p  )  =  temperature  at  another  point  displaced  a  distance  p 
from  x ; 
p  =  autocorrelation  distance; 
<p>  =  a  =  "patch  size". 
Whereas,  the  exponential  function  fitted  the  data  points  and  decayed 
rapidly  as  p  increased  (Figure  2),  it  did  not  have  a  zero  slope  at  the 
origin,  thus,  the  exponential  function  was  not  appropriate  for  very 
small  correlation  distance  (e.g.,  a  few  cm  or  less)  because  of  heat 
conductivity  in  sea  water.   He  proposed  the  Gaussian  function  to 
correct  this  difficulty,  reasoning  that  thermal  conductivity  would 
prevent  any  discontinuity  in  the  temperature  field  [Liebermann,  1951]. 

Liebermann  assumed  that  the  thermal  microstructure  represented  a 
cloud  of  spheres  with  a  Rayleigh  distribution  having  a  mean  size  value 
of  the  correlation  distance  and  therefore  yielding 

R  (p)  =  R  (a)  =  e'1.  (4) 

Using  this  definition,  the  mean  size  of  the  inhomogeneities  ,  a,  was 
found  to  equal  60  cm  in  the  ocean  (Figure  2). 

Piip  [1961]  investigated  the  sound  velocity  microstructure  in  the 
sound  channel  off  Bermuda  with  a  NBS  velocimeter.   It  was  noted  that  a 
number  of  thin  steps  remained  over  periods  of  hours,  some  oscillating 
30  m  in  depth  in  a  three-hour  period.   At  the  same  station  and  same 
depth,  sound  velocity  changes  of  one  m/sec  in  a  few  hours  were  common 
and  "nearly  the  rule".   The  sound  velocity  axis  was  observed  to  change 
150  m  in  depth  and  3.5  m/sec  over  a  three-month  period.   Pipp  showed 
that  the  base  of  the  main  thermocline  (about  1000  to  1300  m)  exhibited 
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relatively  large  and  rapid  changes  in  sound  velocity  in  an  inhomogeneous 
region  in  a  constant  state  of  fluctuation. 
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Figure  2.   Temperature  Fluctuations  and  Associated  Correlation  Function 

(After  Liebermann) 

Small-scale  temperature  and  salinity  structure  was  measured  with  a 
salinity,  temperature,  and  depth  recorder  (STD)  by  Stommel  and  Federov 
[1967]  near  Timor  and  Mindinao  during  July  of  1965.   Their  measurements 
revealed  that  the  vertical  microstructure  actually  consisted  of  very 
thin  high-gradient  layers.   The  thin  layers  were  separated  by  thicker 
layers  of  nearly  homogeneous  (isothermal  or  isohaline)  water,  i.e., 
step-like  in  appearance  on  their  STD  traces.   By  following  a  discrete 
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feature  from  one  cast  to  the  next,  homogeneous  layers  were  found  to  be 
continuous  from  two  to  20  km  horizontally  and  to  range  in  thickness  from 
two  to  40  m  (Figure  3).   A  boot-like  feature  was  found  at  the  bottom  of 
the  mixed  layer  (Figure  4).   At  a  depth  of  200  to  300  m,  the  mean  thermal 
gradient  was  about  0.040C/m.   In  the  thin  steps,  the  gradient  was  less 
than  0.003C/m,  and  in  the  homogeneous  steps,  the  thermal  gradient  was 
greater  than  0.350C/m.   Stommel  and  Federov  concluded  that  the  ocean 
consisted  of  nearly  homogeneous  layers  separated  by  steps  with  extremely 
high  thermal  gradients. 

Tate  and  Howe  [1968],  in  the  summer  of  1966  in  the  Northwestern 
Atlantic,  used  an  STD  to  look  at  the  thermocline  stratification  in  the 
Mediterranean  Water  intrusion  at  a  depth  of  1280  to  1500  m.   They  found 
layering  with  steps  varying  on  the  order  of  0.17  to  0.35C  in  tempera- 
ture, from  0.02  to  0.055  o/oo  in  salinity,  and  from  thickness  of  15  to 
30  m  (Figure  5).   A  step's  average  temperature  change  of  0.25C  produced 
a  a     change  of  -0.044,  while  the  average  salinity  step  of  0.044  o/oo 
produced  a  a  change  of  +0.039.   The  resulting  change  of  a   (-0.005) 
gave  a  slightly  stable  layer  system  in  a  situation  which  was  clearly 
one  of  opposing  salinity  and  temperature  density  gradients. 

Cooper  and  Stommel  [1968],  using  an  STD,  found  the  main  thermocline 
in  the  Sargasso  Sea  to  consist  of  regular  salinity  and  temperature  steps 
of  homogeneous  layers  three  to  five  meters  thick,  separated  by  trans- 
itional layers  10  to  15  meters  thick,  with  temperature  changes  of  0.30 
to  0.50C  and  salinity  changes  of  0.04  to  0.10  o/oo  (Figure  6).   The 
horizontal  extent  of  these  steps  was  between  400  to  1000  meters.   The 
main  thermocline  contained  up  to  100  steps  on  each  cast. 
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Figure  3. a.   Photographic  copy  of  actual  STD  trace  (not  retouched)  of 

station  21  off  Mindinao.   Depth  scale  is  500  m,  each  small 
division  being  5  m.   The  salinity  curve,  labeled  S,  has 
the  smallest  division  being  0.05  o/oo.   The  temperature, 
labeled  T,  has  a  smallest  division  of  0.05C.   The  lamina 
B  is  so  labeled.   Because  the  pens  cannot  both  traverse 
the  same  path  mechanically,  the  temperature  curve  is 
offset  downward  by  5  m  and  always  reads  5  m  too  deep. 
The  salinity  scale  is  not  offset  in  depth.   Uniform 
corrections  for  salinity  of  -0.9  o/oo  and  temperature  of 
0.07C  have  not  been  applied. 

Figure  3.b.   A  composite  chart-graph  for  lamina  B  off  Mindinao.   Each 
square  is  0.05C  wide  (8.5  to  9.0C)  and  50  m  high.   The 
small  arrow  indicates  a  depth  of  400  m  in  each  square. 
Numbers  in  the  lower  left  hand  corners  are  the  station 
numbers.   The  distances  from  stations  18  to  25  and  26  to 
35  are  about  10  nm.   For  lamina  B,  the  large  geographical 
extent  of  this  small  lamina  is  self-evident. 


[Stommel  and  Federov,  1967] 
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Figure  4.   From  the  soundings  in  the  mixed  layer  off  Timor,  gross  dif- 
ferences in  the  properties  in  different  locations  are  self- 
evident.   There  is  considerable  variation  in  vertical 
structure  but  each  sounding  shows  a  nearly  homogeneous 
region  near  the  surface  (station  13)  but  some  stations 
(like  7)  have  shallow  irregular  layers  overlying  the  mixed 
layer. 

[Stommel  and  Federov,  1967] 
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Figure  5 


A  typical  step  in  Mediterranean  Water  [Tait  and  Howe,  1968] 
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Woods  [1968,  1968a]  investigated  the  microstructure  in  the  summer 
thermocline  off  Malta  using  a  temperature-gradient  meter  for  obtaining 
the  thermal  gradient  between  two  thermistors  set  50,  25  or  10  cm  apart 
on  a  vertical  staff.   The  resolution  of  the  thermal  gradient  was  within 
0.01C  per  thermistor  separation  distance.   An  additional  thermistor 
recorded  temperature.   Woods  plotted  profiles  of  the  thermal  gradient 
and  temperature  versus  depth  (Figures  7,  8  and  9).   Using  time- lapse 
photography,  Woods  obtained  pictures  of  streaks  left  in  the  wakes  of 
free-dropping  dye  pellets  for  shear  measurements.   A  dye-packet  array 
tied  between  moored,  submerged  floats  injecting  dye  into  selected  levels 
of  the  thermocline,  indicated  that  the  thermocline  was  divided  into 
layers  a  few  meters  thick.   These  layers  were  characterized  by  weak 
thermal  gradients,  less  than  O.OOlC/m,  and  by  turbulence  with  a  mean 
velocity  of  less  than  one  mm/sec.   The  layers  were  separated  by  sheets 
only  centimeters  thick.   (Sometimes,  he  detected  sheets  up  to  a  meter 
thick  but  they  appeared  as  aggragates  of  the  10  cm  thick  sheets.)   The 
sheets  were  characterized  by  strong  thermal  gradients  up  to  0.05C/cm 
and  little  or  no  turbulence.   Mature  sheets  contained  thin  laminar-flow 
zones  in  their  centers.   Neighboring  layers  had  dissimilar  properties 
indicating  different  salinities,  velocities,  and  turbidities.   This 
suggested  different  trajectories  and  a  very  slow  exchange  of  properties 
across  the  intervening  sheets.   If  turbulence  was  small,  flux  exchange 
was  small;  therefore,  the  exchange  had  to  rely  on  molecular  processes, 
which  were  very  slow  indeed. 

Individual  layers  were  followed  horizontally  for  distances  up  to 
72  km  by  soundings  from  high-speed  launches.   Layers  were  identified  on 
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Figure  8.   Temperature  and  Temperature-Gradient  Soundings  for  Maltese 
Waters,  14  September  1967  [Woods,  1968a]. 

Temperature  0918-0929  local  time,  (i.e.,  GMT  +  1). 

Temperature  gradient  1132-1145  local  time. 
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successive  soundings  for  many  hours,  thus  confirming  the  results  of 
Stommel  and  Federov. 

The  dye  experiments  showed  that  in  a  thermocline  sheet  there  was  a 
central  laminar-flow  region  only  a  few  centimeters  thick.   Dye  was  spread 
in  this  region  by  internal  waves  to  form  a  "broad  carpet  of  color".   This 
colored  sheet  made  it  possible  to  photograph  the  internal  waves  as  they 
crossed  it.   Wave-lengths  between  5  and  10  m  were  dominant  and  appeared 
in  long,  coherent  trains,  each  associated  with  a  single  sheet.   They  had 
heights  up  to  1.0  m,  a  phase  speed  of  2.5  cm/sec  and  periods  of  about 
5  min.   These  waves  only  affected  a  single  sheet,  being  much  shorter  than 
the  long  internal  waves  which  move  the  entire  thermocline. 

A  sheet  was  defined  as  stable  when  it  appeared  smooth  everywhere; 
this  indicated  laminar  flow.   When  a  wave  passed  over  a  sheet,  shear  had 
its  greatest  effects  at  the  crest  and  trough.   Wave  shear  and  drift 
shear  interacted  either  cancelling  or  reinforcing  each  other  causing  a 
sheet  to  become  unstable  and  patchy. 

Wavelets,  about  75  cm  in  length,  caused  banding  to  appear  parallel 
to  the  crests  and  troughs.   These  wavelets  generally  formed,  grew,  and 
broke  in  a  period  of  less  than  two  minutes.   Upon  breaking  in  classical 
fashion,  a  second  smaller  breaker  was  generally  thrown  forward.   The 
wavelets  had  a  maximum  height  of  about  20  cm.   After  breaking,  a  turbu- 
lence region  remained  for  about  5  minutes,  and  the  scar  in  the  sheet 
remained  for  several  hours.   These  scars  were  a  common  feature  in 
thermocline  sheets.   The  thinnest  sheet  observed  to  become  unstable  had 
a  thickness  of  3  cm;  the  thickest  was  30  cm.   The  great  majority  were 
between  8  and  15  cm. 
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Transient  the rmoc lines  [Tully,  1964]  which  occur  sporadically  between 
the  surface  and  the  top  of  the  thermocline  were  considered  as  further 
examples  of  sheets.   Woods  divides  the  thermocline  into  a  half-dozen  or 
so  layers  of  low  shear  and  moderate  thermal  gradient  separated  by  thin 
liminar  flow  sheets  on  intense  shear  and  thermal  gradient.   He  discussed 
the  important  aspect  of  heat  transfer  through  such  a  layered  ocean.   He 
suggested  that  the  molecular  heat  transfer  through  the  high  gradient 
stable  layers  was  negligible,  but  that  heat  transfer  does  take  place 
through  the  turbulent  scars  from  one  layer  to  another  as  an  effect  of  the 
breaking  of  the  internal  waves. 

The  use  of  the  Ramsey  SVTP  (sound  velocity,  temperature,  and  pressure) 
probe  provides  a  digital  output  of  these  parameters.   This  instrument  has 
an  accuracy  of  +  0.01C  and  +  0.25  percent  of  depth  with  a  precision  of 
+  0.08  m/sec  in  sound  velocity  [Lovett,  1968].   In  view  of  the  plus-or- 
minus-one  count  for  electronic  counters,  small-scale  features,  i.e., 
less  than  +  0.01C  and  +  1  m  in  depth,  are  ambiguous.   Lovett  [1968] 
discussed  the  thermal  gradients  determined  with  this  instrument  (Figure 
10).   The  10  m  depth  increments  he  used  eliminated  some  microstructure ; 
however,  using  5  m  increments  introduced  random  readout  errors. 

Ingham  [1968]  used  an  STD  probe  to  investigate  the  mixed  layer  in 
tropical  waters.   The  vertical  salinity  and  temperature  structures 
forming  the  "isothermal"  and  "isohaline"  layers  were  generally  analogous 
in  the  near-surface  layers.   The  difference  in  the  depths  of  these  two 
layers  exceeded  7.5  m  in  45  percent  of  the  casts  and  17.5  m  in  31  percent 
as  shown  by  Figure  11  [Ingham,  1968].   Microstructure  effects  on  sound 
velocity  in  the  near-surface  layers  would  depend  on  salinity  in  some  ocean 
regions . 
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Figure  10.   Temperature  Gradient  per  10  m,  Mean  Gradient  per  100  m,  and 
Average  Deviation  from  Mean  versus  Depth  from  Ramsey  Probe 
Data  [Lovett,  1968]. 
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Figure  11.   Histogram  of  Depth  Differences  AZ  of  the  Isothermal  and 
Isohaline  Layer  Depths  [Ingham,  1968]. 
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Pingree  [1969]  using  an  STD,  recorded  analogue  signals  on  tape.   He 
made  fine-scale  plots  of  salinity  and  temperature.   In  particular,  he 
investigated  the  stability  of  the  deep  water,  i.e.,  about  500  to  2000  m 
(Figure  12). 

Osborn  [1969]  using  a  free-falling  temperature  gradient  meter 
obtained  excellent  profiles  of  temperature  gradient  (Figure  13).   Off 
San  Diego,  the  horizontal  extent  of  the  layers  was  greater  than  750 
meters  in  the  seasonal  thermocline,  but  only  a  few  hundred  meters  at 
depths  greater  than  400  meters.   He  noted  the  layers  were,  in  fact, 
made  up  of  even  finer  layers  with  a  thickness  of  only  a  few  centimeters. 
His  average  sheet  was  about  one  meter  thick. 

Recent  work  by  Neshbya,  Neal,  and  Denner  [1969]  using  the  expendable 
bathythermograph  manufactured  by  General  Motors  Defense  Laboratories, 
obtained  records  of  temperature  steps  (Figures  14,  15  and  16).   While 
on  Ice  Station  T-3  in  the  Arctic  (near  84-38N,  128-22W) ,  they  lowered 
XBT  probes  by  hand.   The  Ice  Station  results  are  the  first  XBT  micro- 
structure  records  without  ship  motion. 

B.   APPLICATIONS 

Urick  and  Searfoss  [1949]  calculated  theoretically  that  a  sound  ray 
would  be  bent  through  an  angle  of  only  3.7  sec  by  the  average  inhomo- 
geneity  in  the  mixed  layer  but  as  much  as  2 . 1  min  of  arc  in  the 
thermocline . 

Sheehy  [1950]  measured  sonar  signal  fluctuation  as  functions  of  depth 
and  range.   He  found  the  coefficient  of  variation,  defined  as  the  stand- 
ard deviation  of  the  amplitudes  expressed  as  a  percentage  of  the  mean 
amplitude,  varied  as  the  square  root  of  range  and  did  not  appear  to  be  a 
function  of  depth. 
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Figure  12.   T-S  Curve  from  1500  to  1800  m  with  a  Value  Every  Meter, 
Circles  Represent  Increments  of  2  m  in  Depth  [Pingree, 
1969]. 
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Pigure    13.      T«»p«r«ture  Gr*di«nt-D«pth  V«lu««    [Ofkora,    1949] . 
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Figure  14.   Vertical  Profile  of  Temperature  Under  Ice  Station  T-3 
(84-38N,  128-22W)  19  March  1969.   [Neshbya,  Neal  and 
Denner,  1969]. 
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Figure  15.   Time  Series  from  2145Z  17  March  to  2033Z  18  March  1969 
[Neshbya,  Neal  and  Dermer ,  1969]. 
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Figure  16.  Time  Series  with  High  Amplification  Taken  from  T-3  in  March 
1969.  Temperature  increases  from  right  to  left.  [Neshbya, 
Neal  and  Denner,  1969]. 
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Liebermann  [1951]  described  the  influence  of  the  microstructure  on 
acoustic  propagation   He  wrote  the  sharpness  of  temperature  boundaries 
determines  influence  of  the  microstructure  on  acoustic  propagation. 
When  temperature  changes  occur  within  a  distance  of  less  than  one 
acoustic  wavelength,  the  temperature  boundaries  may  reflect  sound 
energy  depending  upon  the  size  and  shape  of  the  "patch"   The  tempera- 
ture microstructure  may  reflect,  refract,  scatter,  or  focus  acoustic 
waves,   The  microthermal  structure  was  considered  to  be  made  up  of 
randomly  sized  "spheres",  spaced  randomly  throughout  the  medium,  with 
a  scattering  cross  section  dependent  upon  the  autocorrelation  of  the 
temperature  inhomogeneit ies ,  R  (p) .   Comparing  derived  scattering  cross 
section  with  actual  scattering  in  the  ocean  (i.e.,  volume  reverberation) 
Liebermann  concluded  the  observed  reverberation  was  higher  than  one 
would  expect  from  thermal  inhomogeneit ies   He  attributed  volume 
reverberation  to  biological  scatters. 

Liebermann  found  that  refraction  effects  of  the  inhomogeneit ies 
resulted  in  warping  of  the  wave  fronts,  and  caused  the  wave  front  to 
be  a  complex  surface  exhibiting  both  convergence  and  divergence.   The 
temporal  motion  of  the  inhomogeneit ies  resulted  in  sound  "scintillation" 
at  the  receiver. 

Mintzer  [1953]  found  Sheehy's  coefficient  of  variation  to  be  depend- 
ent upon  the  three-halves  power  of  the  range  rather  than  the  one-half 
power  of  range . 

Skudrzyk  [1957]  discussed  the  focusing  and  defocusing  effect  within 

2 
the  focusing  range  rQ  given  by  r  ■  ka  where  a  is  the  "patch"  radius 

and  k  is  the  acoustic  wavenumber   At  ranges  greater  than  r  ,  in  the 

interference  range,  phase  interference  results  (Figures  17  and  18). 
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Figure  17.   Illustration  of  the  Various  Types  of  Acoustic  Scattering 
Due  to  a  Thermal  Patch  in  the  Sea  [Whitemarsh,  Skudrzyk 
and  Urick,  1957]. 


Figure  18.   The  Phase  Angle  Between  the  Incident  and  the  Scattered 
Pressure  and  the  Physical  Nature  of  the  Scattering  as  a 
Function  of  the  Range  [Skudrzyk,  1957]. 
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Using  this  concept,  Skudrzyk  (Figure  19)  shows  the  focusing  range  r  as 
a  function  of  frequency  and  patch  size.   (For  example,  a  50  kHz  frequency 
sonar  has  a  focusing  range  of  75  meters,  if  the  patch  size  is  assumed  to 
be  60  centimeters.)   For  finding  patch  size,  Skrudzyk  used  the  Urick- 
Searfoss  [1948]  data  to  show  patch  size  was  twice  the  depth  (Figure  20). 
Further,  Skudrzyk  started  theoretical  investigation  with  spheroidal 
patches,  which  he  found  to  have  a  considerably  stronger  focusing  effect 
than  the  spherical  patches. 

Whitemarsh,  Skudrzyk,  and  Urick  [1957]  did  additional  work  with  the 
spheriodal  theory.   The  correlation  functions  previously  used  were  dis- 
cussed.  It  was  concluded  that  the  present  correlation  functions  were 
not  suitable  and  either  had  to  be  replaced  or  redefined.   The  authors 
argued  if  the  temperature  microstructure  was  produced  by  turbulent 
motion  in  the  sea,  then  the  Kolmogorov  Law  for  statistical  equilibrium 
must  hold.   This,  they  show,  was  in  excellent  agreement  with  their  data. 
Lending  support  to  continued  studies  such  as  are  presented  in  this 
thesis,  their  final  conclusion  states  that  the 

"...  greatest  gap  seems  to  lie  in  the  microthermal 
statistics  of  the  ocean.   If  these  were  known  better 
than  they  are  now,  apparently  the  amplitude  statistics 
of  the  fluctuations  of  sound  transmitted  through  the 
sea  could  be  predicted  roughly,  but  with  some  assurance."1 

Works  by  Peterson  [1965],  Barakus  [1968]  and  DiNapoli  [1969]  have 

considered  the  normal  mode  approach  to  the  interference  range  solution 

for  the  temperature  inhomogeneities  and  have  also  investigated  the 

effects  of  internal  waves  on  sound  in  stratified  mediums. 


%/hitemarsh,  D.  C,  Skudrzyk,  E.,  and  Urick,  R.  J.,  "Forward  Scattering 
of  Sound  in  the  Sea  and  Its  Correlation  with  the  Temperature  Micro- 
structure,"  J.  Acous.  Soc.  Am. .  v.  29,  p.  1141,  October  1957. 
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Figure  19.   Focusing  Range  (meters)  as  a  Function  of  Frequency  and 
Patch  Size  (meters)  [Skudrzky,  1957]. 
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Figure  20.   Patch  Diameter  Versus  Depth  [after  Skudrzyk,  1957] 
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C.   INADEQUACIES  OF  THE  PRESENT  THEORY 

Since  the  thermal  microstructure  is  not  related  simply  to  the  gross 
thermal  features  measured  on  MBTs  ,  there  is  presently  no  easy  way  to 
predict,  on  a  day-to-day  basis,  the  effects  of  microthermal  structure 
on  the  propagation  of  sound  in  the  sea.   The  theories,  so  far,  have 
developed  in  two  classes. 

The  first  theories  [Liebermann,  1951;  Mintzer,  1953,  1953a,  1954; 
and  Skudrzyk,  1957]  concerned  spherical  and,  more  recently,  spheriodal 
shaped  inhomogeneities  throughout  the  ocean.   Part  A  above  has  dis- 
cussed the  present  knowledge  of  microstructure  showing  that  the  long 
thin  sheets  and  layers  do  exist  as  microthermal  structure.   These  first 
acoustical  theories  do  not  apply  to  what  is  now  described  as  the  micro- 
structure  of  the  ocean. 

The  theoretical  treatment  by  Whitemarsh,  Skudrzyk  and  Urick  [1957], 
considering  speriodal  patches  and  Kolomogorov' s  turbulence  theories, 
could  not  find  any  real  theoretical  explanation  of  what  happens  in  the 
sea  when  sound  passes  through  an  inhomogeneity .   The  investigations 
using  normal  modes  have  not  produced  any  practical  operational  results. 

Lee  [1961]  using  ray  path  theory,  demonstrated  the  effect  of  a 
typical  sinusiodal  internal  wave  on  sound  intensity  in  a  three-layered 
model.   For  the  medium  with  an  internal  wave,  intensity  variations  up 
to  22  db  were  found  in  distances  of  less  than  one  internal  wavelength, 
as  shown  by  Figures  21  and  22.   In  contrast,  the  same  medium,  without 
an  internal  wave,  produced  a  5  db  variation,  including  a  2.5  db 
spreading  loss,  over  the  same  distance. 

Another  study  investigation  the  effects  of  sound  velocity  micro- 
structure  was  conducted  by  FNWC  [Wolff,  Tatro  and  Megehee ,  1967].   A 


35 


•itrtMCi  in  rur 


Figure  21.   Diagram  of  Sound  Rays  through  a  Medium  Which  has  an 
Internal  Wave  on  the  Thermocline  [Lee,  1961]. 
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Figure  22.   Sound  Level  in  the  Medium  with  an  Internal  Wave  on  the 

Thermocline.  The  db  reference  level  is  that  corresponding 
to  a  sound  level  of  60  db  at  1  foot  from  the  directional 
source  (X  equals  0,  Z  equals  10  ft)  along  the  horizontal.  ■ 
The  field  is  contoured  at  intervals  of  2.5  db  [Lee,  1961]. 
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sophisticated  ray-tracing  program  [Ayers ,  Wolff,  Carstensen  and  Ayers ,  1966] 
was  used  for  the  analysis.   Three  complete  ray  trace  runs  were  made  to 
the  first  convergence  zone  using  120,  54,  and  20  digitized  sound  velo- 
city points  from  one  SVTP  trace  containing  sound  velocity  microstructure . 
The  120-point  sound  velocity  profile  preserved  all  the  "wiggles"  or 
microstructure,  the  54-point  profile  used  values  at  every  50  meter  of 
depth,  and  the  20-point  profile  used  the  standard  oceanographic  depth 
values  for  sound  velocity  (Figure  23).   The  range  to  the  first  conver- 
gence zone  agreed  within  150  yards  for  all  three  cases.   Based  on  the 
ray  trace  runs  only,  it  was  concluded  the  sound  velocity  microstructure 
was  "essentially  self-cancelling"  and  did  not  affect  the  convergence 
zone  parameters  [Wolff,  Tatro,  and  Meghee ,  1967].   Propagation  loss  and 
reverberation  effects  were  not  reported. 
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Figure  23 


A  Portion  of  the  Three  Different  Sound  Velocity  Profiles 
Used  by  FNWC  for  Ray  Tracing   [after  Wolff,  et  al.]. 
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III.   MEASURING  MICROTHERMAL  STRUCTURE 

A.   THE  BATHYTHERMOGRAPH 

1.  Mechanical  Bathythermograph 

At  the  present  time  temperature  is  the  most  readily  observed 
variable  in  the  upper  few  hundred  meters  of  the  ocean.   A  major  con- 
tributor to  this  situation  is  a  singularly  successful  oceanographic 
instrument  known  as  the  mechanical  bathythermograph  or  MBT.   Subsequent 
to  its  introduction  by  Rossby  and  Montgomery  [1934]  and  perfection  by 
Spilhaus  [1938],  the  MBT  has  become  one  of  the  basic  instruments  of 
almost  all  oceanographic  agencies  because  of  its  simplicity  of 
operation,  reliability  and  relatively  low  cost.   The  result  has  been 
that  a  large  number  of  bathythermograms  have  been  collected. 

Researchers  during  World  War  II,  using  the  MBT,  observed  micro- 
thermal  structure  on  their  traces.  No  theory  had  been  formulated  which 
took  into  account  the  effects  caused  by  microthermal  structure  on  under- 
water sound  transmission.  It  was  recognized,  however,  that  temperature 
microstructure  must  be  considered  in  any  explanation  of  observed  under- 
water sound  transmission  [National  Defense  Research  Committee,  1946]. 

2 .  The  Expendable  Bathythermograph 
a.   Capabilities 

Presently  the  U.  S.  Navy  is  using  the  XBT  system  manufactured 
by  the  Sippican  Corporation  of  Marion,  Massachusetts.   This  system  con- 
sists of  an  expendable  shipboard  cannister  and  probe,  shipboard  launcher, 
and  a  strip  chart  recorder.   Data  so  obtained  are  submitted  to  FNWC  and 
were  used  in  this  thesis. 
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The  original  Bureau  of  Ships'  specifications  [Denner,  1966] 
for  an  XBT  system  are  given  in  Table  I. 


TABLE  I 


Temperature  range 
Temperature  accuracy 
Depth  range 
Depth  accuracy 

Ship's  speed  when  launching 


28F  to  96F 

+  0.4F 

0  to  1000  ft 

+  2$  or  15  ft,  whichever 

is  greater 
+  63$  of  a  temperature 

step  after  3  ft  depth 

change 


The  expendable  Sippican  probe  currently  used  has  a 
thermistor  sensor  in  its  nose.   The  probe  is  connected  electrically 
to  the  recorder  by  two  conductor  wires  which  pay  out  from  both  the 
probe  and  the  canister.   The  cannister  in  the  launcher  is  connected 
to  the  recorder  by  a  cable.   There  are  no  splices  in  the  wire.   The 
probe's  weight  is  a  critical  factor  in  determining  its  fall-rate. 
The  Sippican  system  was  found  to  be  within  +  3  gms  both  with  and 
without  its  internal  copper  wire.   (The  deviation  was  from  0.3$  to 
0.64$  of  the  total  weight  of  the  probe.)   The  observed  average  fall- 
velocity  of  the  probe  was  20.3  +  0.2  ft/sec.   A  probe,  with  its 
internal  wire,  had  an  average  terminal  velocity  of  19.9  +  0.7  ft/sec. 
Terminal  velocity  was  reached  at  18.0+0.5  ft  [Gouzie,  Sanders  and 
Littlehale,  1966]. 

The  present  XBT  strip  chart  recorder  is  a  Mark  2A  model 
using  pressure  sensitive  paper  for  recording  the  thermal  profile.   The 
paper  moves  vertically  on  the  recorder  at  a  rate  of  3.25  in/min.   In  one 
minute,  the  XBT,  falling  at  20  ft/sec,  will  have  reached  a  depth  of  1200 
ft.   Thus,  1200  ft  will  be  recorded  on  3.25  inches  of  chart  paper,  or 
0.0027  of  an  inch  of  chart  paper  is  used  for  every  foot  of  fall. 
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The  chart  paper  records  horizontally  from  28F  to  96F  in 
7.0  inches,  or  about  0.103  of  an  inch  of  chart  paper  per  1.0F.   Table 
II  shows  the  resolution  capabilities  of  the  paper  using  different 
recording  intervals. 


TABLE   II 

Resolut] 

Lon 

Chart   Spacing    (inches) 

Temperature    (F) 

Depth    (ft) 

0.01 

0.0971 

3.690 

0.005 

0.0476 

1.045 

0.001 

0.0097 

.369 

The  response  time  for  the  thermistor  is  110  milliseconds 
nominal  and  130  milliseconds  maximum  [Choate,  1970  and  Demeo,  1968]. 
A  time  constant-temperature  change  plot  shows  that  98.2  percent  of  a 
temperature  change  has  been  recorded  in  10.4  feet  of  fall,  83  percent 
recorded  in  5 . 2  feet  of  fall,  68  percent  in  two  feet  of  fall  and  35 
percent  in  one  foot  of  fall  (Figure  24). 

The  thermistor,  because  of  its  finite  size,  has  thermal 
inertia  which  causes  it  to  act  as  a  low  pass  filter.   As  a  result, 
fine  structure  is  lost  or  at  best  smoothed  out.   The  original  trace 
could  be  reconstructed  if  the  filtering  characteristics  (of  both  phase 
and  amplitude)  of  the  thermistor  were  known.   These  could  be  obtained 
experimentally . 

b.   Data  Handling 

Present  procedures  for  handling  XBT  data  are  indicated  by 
the  log  sheets  which  accompany  the  XBT's  in  the  shipping  container. 
First,  a  BATHY  or  BAXBT  message  is  required  to  be  submitted  to  the 
appropriate  addressee  using  a  PRIORITY  precedence  message  [National 
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Oceanographic  Data  Center,  1968,  1968a].   Upon  receipt  of  the  BATHY  or 
BAXBT  message  the  receiving  agency  uses  this  as  synoptic  data  for  sub- 
jective ocean  forecasts  in  his  locality.   The  BATHY  or  BAXBT  is  for- 
warded through  the  Naval  Environmental  Data  Network  to  FNWC  in  Monterey. 
FNWC  uses  the  message  reports  to  update  its  ocean  fields  for  objective 
forecasting  in  various  oceanographic  and  meteorological  products. 

Upon  completion  of  a  cruise  the  ship  which  took  the  XBT 
cast  mails  the  original  trace  and  log  sheet  to  FNWC  for  digitizing. 
The  original  XBT  trace  is  digitized  for  use  by  FNWC,  and,  when 
requested,  the  results  of  the  digitizing  are  provided  to  the  XBT 
submitting  activity.   Upon  completion  of  digitizing,  the  original  XBT 
traces  are  forwarded  to  the  National  Oceanographic  Data  Center  (NODC) 
for  storage  in  their  archives.   NODC  also  receives  a  copy  of  the  digital 
data  on  tape. 

Present  procedures  at  FNWC  use  DIGITBT  and  several  other 
specialized  XBT  data  handling  programs.   New  XBT  digitized  data  are 
stacked  on  the  XBT  Master  File  (presently  at  a  rate  of  about  500  traces 
per  week) . 

B.   OTHER  SENSORS 

Sensors  for  measuring  microthermal  structure,  other  than  those 
discussed  elsewhere  in  this  thesis,  are  numerous. 

From  a  fixed  platform,  like  the  NEL  Oceanographic  Tower,  thermistor 
strings  and  arrays  are  used.   Bottom  devices,  such  as  thermistors, 
attached  to  the  shore  or  a  tower  by  a  wire,  have  been  employed  to  give 
continuous  temperature  records. 
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From  submarines,  thermopiles,  thermocouples,  thermistors,  platinum 
resistance  thermometers,  and  thin  film  thermometers  have  been  employed. 
Almost  any  oceanographic  electronic  or  mechanical  thermometer  could  be 
used  on  a  submarine. 

C.   DIGITIZING  PROCEDURES 

At  the  present  time,  all  XBT  traces  are  submitted  to  Fleet  Numerical 
Weather  Central  (FNWC)  in  Monterey,  California,  the  only  digitizing 
center  for  these  strip  chart  analogue  traces  [National  Oceanographic 
Data  Center,  1968,  1968a;  and  Bauer,  1969]. 

At  the  Point  Pinos  Annex  of  FNWC,  two  Ca lma  Model  480  flat -bed 
digital  recorders  are  used  to  digitize  XBT  traces  in  0.01  inch 
increments.   These  digital  recorders  are  operated  in  a  mode  which  makes 
the  recorded  intervals  dependent  on  the  tracing  speed.   The  Calma 
output  recorded  on  digital  magnetic  tape  containing  binary  mode 
characters  at  a  density  of  556  bits  per  inch  (BPI)  with  even  parity 
and  variable  length  records. 

FNWC  program  DIGITBT,  written  by  Lieutenant  N.  L.  Perkins,  USN, 
in  FORTRAN  IV  for  FNWC ' s  CDC  6500  computer,  is  used  currently  for 
the  reduction  of  data  from  XBT  traces.   FNWC  tape  input,  tape  output, 
and  some  other  subroutines  of  DIGITBT  are  written  in  COMPASS  language. 
In  order  to  analyze  microthermal  structure,  DIGITBT  was  extensively 
modified.   The  new  program,  referred  to  as  MICROXBT,  is  included  as 
Appendix  A.   The  basic  procedures  for  DIGITBT  and  MICROXBT  programs 
are  the  same . 

DIGITBT  and  MICROXBT  programs  both  convert  the  X  and  Y  measure- 
ments of  the  digitizer  into  depth  (X)  versus  temperature  (Y)  values. 
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Corrections  are  made  for  nonlinearities  produced  by  the  XBT  recorder  in 
both  the  depth  and  temperature  directions.   All  shipboard  XBT  recorders 
are  assumed  to  produce  error-free  traces.   The  62 . OF  temperature  line 
drawn  on  the  trace  when  the  probe  is  loaded  is  used  as  a  baseline  for 
digitizing.   No  correction  is  made  if  this  baseline  temperature  line 
is  not  the  same  as  the  62. OF  line  printed  on  the  trace  paper.   Ah 
eye-interpolated  value  of  the  XBT-drawn  baseline  temperature  on  the 
printed  temperature  grid  is  noted  as  a  part  of  the  identification 
record. 

The  rigid  procedures  for  digitizing  require  that  exactly  one 
identification  (ID)  record  and  two  digitized  trace  records  for  each 
XBT  be  recorded  on  the  tape.   The  ID  record  information  contains  the 
data  (day,  month,  year);  time  (hour  and  minute);  latitude  and  longi- 
tude (degree,  minute,  hemisphere);  baseline  temperature  (to  nearest 
0.1C  or  0 . IF  depending  on  the  grid  type);  printed-grid  type;  and  maxi- 
mum depth  for  the  probe  (coded  as  15  for  1500  ft  XBT  or  25  for  the 
2500  ft  XBT).   Each  trace  record  contains  pairs  of  increments  (AX,  AY) 
read  from  the  baseline  to  the  true  surface  temperature,  a  point  "flag" 
and  then  pairs  of  values  (AX,  AY)  along  the  trace  to  a  cut-off  depth. 
Each  trace  is  digitized  twice.   The  complete  procedure  is  described 
by  Bauer  [1969] . 

To  process  digitized  XBT  data  using  the  DIGITBT  program,  the  two 
trace  records  are  compared  point-by-point  in  a  deviation  test.   If 
the  X  and  Y  paired  values  (i.e.,  depth  and  temperature)  differ  by  more  than 
0.035  inches,  DIGITBT  rejects  the  whole  XBT  record  and  the  XBT  trace  must 
be  completely  re-digitized.   If  the  XBT  record  passes  the  deviation  test, 
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the  two  traces  are  averaged,  point-by-point.   The  resultant  trace  is 
used  for  the  rest  of  the  program. 

In  the  following  the  term  "inflection  point"  will  be  used  to  pre- 
serve consistency  with  the  author  of  DIGITBT.   Inflection  points  are 
those  points  which  meet  the  programs' criteria  as  significant  points 
which  are  the  programs  output. 

Both  DIGITBT  and  MICROXBT  select  the  first  and  last  observations 
as  inflection  points  in  the  output  list.   The  two  programs  differ  at 
the  inflection  point  choice  stage.   Other  observations  in  DIGITBT  are 
screened  in  groups,  and,  for  each  group,  the  single  observation  is 
found  which  is  most  different  from  a  linear  interpolation  between  the 
two  end  points  of  this  group  of  observations.   Each  observation  in  a 
group  is  measured  to  find  how  out-of-line  it  is  in  temperature  and  in 
depth,  and  the  smaller  of  these  two  is  considered  the  out-of-line 
measure.   All  observations  in  the  group  are  then  compared  to  find  the 
one  with  the  largest  out-of-line  measure.   If  that  value  exceeds 
0.0150  for  an  inch,  that  observation  is  selected  for  output  in  DIGITBT; 
otherwise,  the  group  expands  to  include  the  next  value  and  the  inflec- 
tion test  is  repeated  until  an  inflection  point  is  found. 

The  selection  of  inflection  points  was  modified  in  MICROXBT  to  use 
a  minimum  of  0.0001  of  an  inch  for  the  out-of-line  measure.   This  pro- 
duces an  output  containing  every  inflection  point,  which  is  nonlinear 
between  two  consecutive  points  in  the  digitized  record.   Any  observations 
which  cannot  be  reconstructed  by  linear  interpolation  are  not  required 
and  are  removed  from  the  output  listing  in  both  DIGITBT  and  MICROXBT. 

DIGITBT  converts  depth  units  from  digitized  inches  to  ft  of  depth 
by  linear  interpolation  using  a  depth  table  and  from  ft  to  m  by  formula. 
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Temperature  units  are  converted  from  digitized  inches  to  F  by  a  tempera- 
ture table  and  from  F  to  C  by  formula.   These  depth  and  temperature 
tables  were  constructed  to  account  for  the  nonlinearities  between  the 
recorder  and  the  actual  trace.   Values  of  depth  are  rounded  to  the 
nearest  m  and  temperature  to  the  nearest  0.01C  in  DIGITBT.   MICROXBT 
has  eliminated  the  ft  to  m  and  F  to  C  conversions.   Depths  are  rounded 
to  the  nearest  0.1  ft  and  temperature  to  the  nearest  0.001F. 

After  the  conversion  from  ft  to  the  nearest  m,  DIGITBT  removes, 
for  a  second  time,  any  of  the  rounded-off  values  which  are  nonlinear 
between  two  neighboring  points.   If  depth  values  appear  as  duplicates, 
one  value  is  omitted  in  the  output  listing.   In  MICROXBT,  there  are  no 
duplicate  values  because  0.1  ft  in  depth  is  well  beyond  the  hundredth 
of  an  inch  digitizing  capability  of  the  digitizer.   Thus,  by  reading 
depths  in  tenths  of  a  ft  and  temperature  in  thousandths  of  a  degree  F, 
MICROXBT  looses  no  values  from  duplicate  depths,  and  the  possibility  of 
introducing  false  linearities  is  reduced  considerably.   For  these  rea- 
sons, MICROXBT  gives  a  much  finer  digitized  record  than  DIGITBT. 

The  DIGITBT  output  tape  is  written  in  binary  for  the  CDC  6500;  the 
MICROXBT  output  tape  is  written  for  the  IBM  360/67  at  the  Naval  Post- 
graduate School  (NPS).   Appendix  C  contains  the  job  control  cards 
(JCL)  required  for  using  a  CDC  6500  output  tape  from  MICROXBT  with 
external  BCD  at  a  density  of  556  BPI  as  in  input  to  a  separate  micro- 
thermal  analysis  program  at  NPS. 
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IV.   MODELS  FOR  EXAMINING  MICROTHERMAL  STRUCTURE 

All  models  and  methods  used  in  this  thesis  required  the  output  tape 
from  MICROXBT  to  be  used  as  a  direct  input  on  the  seven  channel  tape 
unit.   The  XBT  ID  record  information  was  read  first  and  if  the  recorded 
XBT  was  from  a  desired  geographical  region  the  data  (depth  and  tempera- 
ture) were  analyzed  for  microthermal  structure  at  NPS . 

A  "best"  value  for  the  temperature  structure  was  found  by  using  the 
depth  inflection  points  from  MICROXBT.   For  this  "best"  array,  tempera- 
ture was  linearily  interpolated  to  a  value  for  every  foot  of  depth.   The 
mixed  layer  depth  (MLD)  was  determined,  using  inflection  points  from 
MICROXBT,  as  that  depth  immediately  above  the  depth  where  the  temperature 
was  2 . OF  different  from  that  at  the  surface  [Naval  Weather  Service,  1967] 

Several  microstructure  models  were  tried  in  order  to  perform  an 
objective  analysis  of  the  oceanic  microthermal  structure.   The  use  of 
a  numerical  temperature  derivative  model  was  found  to  give  the  most 
satisfactory  presentation  of  the  microthermal  structure. 

A.   DERIVATIVE  MODEL 

1 .   Mathematical  Model 

The  most  significant  feature  of  the  thermal  structure  is  the 
thermal  gradient  defined  as  the  derivative  of  the  temperature  with 
respect  to  depth  (dT/dZ).   By  plotting  the  gradient  versus  depth  (Z) , 
the  thermal  structure  of  the  water  column  is  emphasized.   Furthermore, 
a  parameter  such  as  the  thermal  gradient  is  related  to  stability, 
entrophy  generation,  turbulent  heat  transport,  and  the  Brunt-Vaisala 
frequency  [Osborn,  1969]. 
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2 .   Experimental  Methods 

There  are  two  methods  for  experimentally  determining  the  thermal 
gradient  presently  in  use. 

a.  Electronic 

A  resistance  thermometer  in  an  electronic  bridge  passing 
through  water  produces  a  signal  which  is  proportional  to  temperature  as 
a  function  of  time.   By  differentiating  this  signal  with  respect  to  time 
(which  is  analogous  to  differentiating  with  respect  to  distance),  the 
values  of  gradient  versus  distance  the  probe  travels  are  found. 

This  method  was  successfully  employed  by  Grant,  Moillet  and 
Vogel  [1968],  in  the  horizontal  direction  for  temperature  and  turbulence 
measurements.   A  thin-film  resistance  thermometer,  mounted  on  the  bow  of 
a  submarine,  measured  temperature  versus  horizontal  distance  traveled. 
The  time  derivative  was  taken  to  produce  a  direct  plot  of  dT/dX,  the 
horizontal  thermal  gradients. 

Osborn  [1969]  used  a  thermistor,  capacitively  coupled  to 
two  differentiating  amplifiers,  to  measure  the  vertical  temperature 
gradient  with  his  free-falling  gradiometers .   His  resulting  output 
(Figure  13)  shows  a  gradient  versus  depth  record  for  the  North  Pacific 
off  the  San  Diego  Trough. 

b.  Two-probe  System 

Temperature  gradients  may  be  found  by  using  two  vertically 
spaced  thermistors .   The  temperature  difference  is  a  measure  of  the  mean 
gradient  between  them.   Such  a  gradient-meter  was  used  by  Woods  [1968, 
1968a]  who  obtained  gradient  at  various  depths  (Figures  7,  8,  and  9). 
This  method  is  limited  to  the  resolving  power  of  the  thermistors  and  is 
a  function  of  thermistor  spacing. 
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3 .   Numerical 

The  numerical  method  is  based  on  the  derivative,  dT/dZ ,  for  the 
gradient  but  uses  the  finite  difference  approximation.   The  derivative 
of  the  temperature  at  any  depth  T(Z.)  with  respect  to  depth  Z  is  defined 


as 


dT  ..  lim   T<V  -  T<Z> 

dZ  '  Z.-Z    Z.  -  Z 

1       1 


(5) 


which  can  be  shown  to  equal: 


dT  _  lim  AT 
dZ  "  AZ-0  AZ 

(6) 

lim  T(Z  +  AZ)  -  T(Z) 

AZ-0  "      AZ 

For  this  thesis,  AZ  =  1  foot,  which  yields  the  finite  difference  approxi- 
mation for  the  vertical  temperature  derivative: 

B.   MISCELLANEOUS  METHODS 

1.   Statistical  Bias  and  Variance 

For  this  method,  a  standard  or  smooth  curve  was  required. 
Several  standard  curves  were  investigated.   The  variance  and  bias  were 
calculated  between  the  standard  and  the  "best"  curves.   The  standard  was 
found  by  computing  the  temperature  values  for  a  series  of  specified  depths 

The  first  standard  curve  investigated  used  MID  and  the  FNWC 
significant  levels  (surface,  100,  200,  300,  400,  600,  800,  and  1200  ft). 
The  temperature  values  were  found  using  the  "best"  temperature  values  at 
those  depths.   Next,  an  array  of  temperature  values  was  found  at  foot 
increments  by  linear  interpolation  between  the  temperature  values  at  the 
significant  levels  or  the  MID.   This  was  designated  the  FNWC  significant 
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level  curve.   Calculations  of  bias  and  variance  did  not  reveal  any  cor- 
relation with  the  microthermal  structure  easily  visible  on  the  XBT  trace. 

The  second  standard  curve  investigated  used  the  temperatures  at 
the  standard  oceanographic  station  depths  (surface,  10,  20,  30,  50,  70, 
100,  150,  200,  250,  300,  400,  500,  and  800  m) .   Since  these  depths  are 
metric,  a  linear  interpolation  was  done  between  the  two  neighboring 
values  on  the  "best"  curve.   Values  of  temperature  for  every  foot  were 
found  by  linear  interpolations  between  the  standard  depth  temperature 
values.   Again,  no  significant  correlation  was  noted. 

The  third  standard  curve  used  thse  standard  depth  temperature 
values  but  used  these  standard  depth  values  as  an  input  to  the  LGINTP 
subroutine  [Denner,  1969].   This  subroutine  uses  the  four  standard  depth 
points  neighboring  the  depth  in  question  to  make  an  "average"  curve  from 
an  upward  and  downward  parabolic  interpolation  and  a  linear  interpolation. 
This  scheme  produces  a  smoothing  of  points  on  the  curve.   This  technique 
is  usually  used  for  interpolating  standard  depth  oceanographic  data 
[Rattray,  1965].   These  values,  after  computing  the  bias  and  variance, 
again  gave  no  correlation  with  the  microthermal  structure. 

The  LGINTP  subroutine  was  not  used  with  the  FNWC  significant 
level  values  for  smoothing  since  the  FNWC  fields  assume  linear  inter- 
polation. 

Other  standards  considered,  but  not  investigated,  included  a 
Gaussian  standard,  using  the  technique  for  approximating  the  temperature 
structure  proposed  by  Grosfils  [1968]  and  an  exponential  standard. 
2 .   Gradient  Distribution  Method 

For  investigating  the  temperature  structure,  the  "best"  tempera- 
ture data  were  separated  into  the  FNWC  significant  levels.   For  each 
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significant  level,  the  temperature  change  from  one  depth  to  the  next 
were  calculated  and  then  sorted  in  numerical  ascending  order  with  the 
IBM  subroutine  SHSORT.   With  these  values  of  temperature  change  per  foot 
or  gradient  (F/ft),  the  temperature  structure  was  now  characterized  by 
the  individual  thermal  gradient  values  and  a  total  temperature  change  AT. 
Because  of  the  range  of  temperature  values  in  the  thermocline,  this  region 
of  the  XBT  trace  was  chosen  for  developing  this  method. 

The  capabilities  of  the  XBT  indicate  a  gradient  interval  of  0.02 
F/ft  is  easily  obtained  with  a  high  degree  of  confidence  [Arthur  D. 
Little,  1966].   Using  0.02  F/ft  for  the  class  interval,  the  frequency 
distribution  of  gradients  was  plotted  for  various  XBT  traces  by  signifi- 
cant levels.   As  expected,  the  greatest  range  in  gradient  values  usually 
occurred  in  the  significant  level  just  below  the  mixed  layer  where  the 
range  of  temperature  values  AT  was  the  greatest. 

This  method  appears  to  give  a  good  indication  of  the  variability 
of  the  thermal  structure  in  the  depth  intervals  investigated.   Several 
examples  of  the  different  types  of  gradient  distributions,  using  histo- 
grams and  a  cumulative  frequency  diagram,  are  shown  in  Figures  25  through 
31.   The  actual  analyzed  XBT  is  also  shown  on  the  left  of  each  figure. 
Looking  at  data  off  San  Diego,  the  first  three  gradient  distribution 
diagrams  (Figures  25,  26  and  27)  show  fairly  similar  cumulative  fre- 
quency distributions.   The  next  two  (Figures  28  and  29)  show  nearly 
identical  distributions  for  water  inside  the  Gulf  Stream.   The  next  two 
gradient  distributions  (Figures  30  and  31)  show  water  in  the  Gulf  Stream's 
frontal  region.   These  last  two  distributions  are  not  identical  for  many 
reasons,  but  mainly  because  the  positive  gradient  region  in  the  second 
example  is  below  the  100  to  200  foot  interval  shown  in  the  gradient 
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distribution  plot.   These  results  indicate  that  the  small  scale  structure 
may  be  water  mass  related.   However,  this  could  not  be  fully  examined  in 
the  time  alloted  for  this  thesis  but  is  recommended  for  further 
investigation. 
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V.   APPLICATION  OF  NUMERICAL  DERIVATIVE  MODEL 

This  chapter  describes  the  application  of  the  numerical  derivative 
model  or  method  (NDM)  to  several  cases  of  real  and  assumed  data. 

A.   BASIC  CASES 

Graphs  are  used  to  show  the  features  of  the  NDM.   For  each  station, 
temperature  versus  depth  (T-Z)  and  temperature  gradient  versus  depth 
(G-Z)  plots  were  drawn  by  a  CalComp  plotter  using  the  computer  program 
with  the  DRAW  subroutine  in  Appendix  B.   The  scales  are  constant  for  all 
plots.   The  grid  squares  were  drawn  for  100-ft  depth  increments  from  the 
surface  to  a  depth  of  900  ft.   On  the  T-Z  plots  (reporduced  digitized  XBT 
data),  temperature  was  scaled  at  10F  per  grid  line.   On  the  G-Z  plots,  the 
gradient  was  scaled  to  O.lF/ft  per  grid  interval,  with  the  axis  offset  so 
that  two-thirds  of  the  plot  was  allocated  to  negative  gradients.   This 
scaling  was  chosen  so  that  each  hundredth  of  an  inch,  the  resolving 
power  of  the  CalComp  plotter,  represented  one  foot  of  depth  and  0.01F 
of  temperature  or  O.OOlF/ft  of  temperature  gradient.   The  plots  included 
here  have  been  reduced  in  size,  but  the  relative  proportions  remain  the 
same  . 

1.   Linear  Gradients 

For  a  temperature  distribution  with  linear  thermal  gradients, 
the  G-Z  plot  is  characterized  by  steps  indicating  a  constant  gradient 
during  a  temperature  change.   For  this  thesis  the  strength  of  a  gradient 
is  defined  as  the  absolute  value  of  the  gradient  (in  F/ft  units)  and  is 
described  as  strong  when  greater  than  0.3F/ft  or  weak  when  less  than 
0.05F/ft.   The  term  signed  strength  is  defined  as  the  true  gradient 
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value  and  will  thus  distinguish  absolute  and  true  gradient  values.   The 
magnitude  of  a  gradient  is  defined  as  the  vertical  extent  of  the  gradient 
measured  in  feet.   The  product  of  magnitude  and  signed  strength  gives  the 
total  temperature  change  for  the  extent  of  the  gradient. 

a.  Layer  Structure 

Figure  32  shows  a  two- layer  temperature  structure  model  with 
a  T-Z  plot  on  the  left  and  a  G-Z  plot  on  the  right.   The  strength  of  the 
upper  layer  is  0.05F/ft  with  a  magnitude  of  200  ft.   The  lower  layer  is 
weaker  (0.02F/ft  compared  to  0.05F/ft)  than  the  upper  step  but  the  magni- 
tude is  greater  (600  ft  to  200  ft). 

For  the  mu It i- layered  model  (Figure  33)  several  linear 
gradient  steps  are  shown.   In  the  gradient  region  between  70  and  90  ft, 
the  gradient  is  strong  and  goes  off  scale  (Figure  33b). 

b.  Step  Structure 

A  step  structure  is  characterized  by  thin  strong  gradient 
regions.   The  model  (Figure  34)  shows  clearly  that  the  four  steps  bet- 
ween 60  and  100  ft  are  stronger  than  any  other  steps  shown.   These  four 
steps  also  have  the  smallest  magnitude  (five  ft)  in  the  figure.   The 
steps  at  200,  300,  and  400  ft  all  have  the  same  strength  but  increase 
in  magnitude  with  depth. 

The  features  of  the  T-Z  trace  (Figure  34a)  are  easily 
described  by  the  G-Z  model.   The  product  of  signed  strength  and  magnitude 
gives  the  total  temperature  change  from  the  surface  temperature  to  a 
given  depth  for  a  step-structure  temperature  profile.   The  temperature 
at  any  depth  could  then  be  found  if  the  sea  surface  temperature  was 

iBj   Thus,  with  magnitude,  signed  strength  and  sea  surface  temperature, 
any  layer,,      Lepped  thermal  profile  could  be  reconstructed. 
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Figure  32.   Two-Layer  Model 
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Plgure  33.   Multilayer  Model 
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Figure  34.   Step  Structure  Model, 
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2 .   Continuous  Curving  Gradients 

Since  finite  values  are  used  for  both  depth  and  temperature,  a 
truly  continuous  curve  cannot  be  drawn.   The  spacing  by  the  CalComp 
plotter  of  every  hundredth  of  an  inch  can  cause  variations  of  this 
size.   Computed  values  of  depth,  temperature  and  temperature  gradient 
were  provided  as  inputs  to  a  CalComp  plotter  at  every  hundredth  of  an 
inch  to  create  a  quasi-continuous  curve,  at  least  to  a  resolution  of  the 
plotter.   It  is  this  quasi-continuous  field  which  accounts  for  the  0.01 
inch  steps  encountered  in  curves  drawn  by  the  plotter  for  any  truly  con- 
tinuous curve.   This  effect  is  also  present  whenever  the  plotter  is 
drawing  any  line  not  directly  in  the  X  or  Y  axis  directions.   Accordingly, 
any  steps  shown  of  this  size  could  be  considered  as  possible  adjustments 
made  by  the  plotter  or  data  to  adjust  to  the  hundredth  of  an  inch  grid. 

Two  examples  of  continuous  gradients  are  given  to  show  the 
resulting  distributions  of  gradient  with  respect  to  depth. 

a.   Gaussian 

An  example  of  a  Gaussian  profile  (Figure  35)  shows  a  normal 
distribution  of  temperature  with  depth  for  the  function,  like 


Z      1   2 
-  i-  x 


20  r       '  2 

T(z)  - 80  -  jr,  }  e       dx 


(8) 


\ 


where:     x  =  a  dummy  variable 

80  =  sea  surface  temperature. 
The  G-Z  curve  highlights  the  changes  in  the  T-Z  curve.   The  increase  in 
gradient  strength  is  very  clear,  reaching  a  maximum  200  feet,  whereas 
the  inflection  point  of  the  T-Z  plot  is  difficult  to  find  by  visual 

ii.   The  summation  of  signed  strength  values  at  every  foot  of 
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Figure  35.  Gaussian  Model. 
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depth  from  the  surface  down  gives  a  total  temperature  change  which  is 
analogous  to  the  product  of  magnitude  and  signed  strength  with  all 
magnitudes  set  equal  to  one.   With  this  summation  and  a  given  sea  surface 
temperature  (SST) ,  a  continuous  temperature  profile  can  be  reconstructed, 
b.   Exponential 

Using  an  exponential  expression  for  the  temperature  distri- 
bution 

T(Z)  =  60  -  20e-nZ/100°  (9) 

where  : 

60  =  SST, 
the  T-Z  (Figure  36a)  and  G-Z  curves  (Figure  36b)  appear  as  compliments 

of  each  other.   The  derivative  of  T(Z)  with  respect  to  Z  gives: 

,n-   „  -nZ/1000 

**S21  =  us.  (l0) 

dZ     50  KW) 

which  is  the  equation  of  the  G-Z  plot. 

3 .  Superpositions 

In  the  example  of  layers  superimposed  on  an  exponential  model 
(Figure  37),  the  conditions  for  the  exponential  (9)  hold  for  the  depth 
intervals:   surface  to  40  ft,  220  to  340  ft  and  460  to  800  ft,  with 
various  layers  at  the  remaining  depths.   Isothermal  layers  exist  below 
the  weak  surface  exponential  gradient  layer  and  also  in  the  140  to  180  ft 
depth  layer.   From  the  T-Z  plot,  the  layers  at  340  and  400  ft  are  not  as 
evident  as  they  are  in  the  G-Z  plot. 

4 .  Invers  ions 

A  model  with  two  inversions  and  several  linear  steps  (Figure  38) 
shows  the  positive  gradient  regions  clearly.   The  second  deeper  inversion 
learly  a  weak  gradient,  but  it  is  still  discernible  in  Figure  38b. 
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Figure  36.   Exponential  Model 
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Figure  37.   Exponential  Model  with  Layers  Superimposed 
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Figure  38.     Multi-Step  Model  With  Two  Inversions 
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With  inversions,  summing  the  signed  strength  values  for  each  foot  gives 
the  temperature  change  from  the  surface  and  adding  this  value  to  SST 
yields  the  actual  temperature  at  any  depth. 

B.   OBSERVATIONAL  CASES 

Several  different  oceanic  regions  were  selected  to  establish  fea- 
tures of  this  model.   The  oceanic  regions  explored  were  along  the  Atlantic 
Coast*  in  the  Andaman  Sea,  Eastern  North  Pacific,  and  the  Gulf  Stream. 

1.  Atlantic  Coastal  Region 

Stations  in  Massachusetts  Bay  occupied  during  December  1969 
were  used.   The  stations  were  within  a  radius  of  13  miles.   Water  depths 
were  between  280  and  420  feet.   Lack  of  strong  gradients  in  all  seven 
traces  shows  that  the  water  was  probably  well  mixed  (Figures  39  through 
42).   A  positively  signed  strength  in  the  upper  300  feet  was  common  in 
the  traces.   Inversions  (Figure  40)  are  quite  clear  in  both  the  T-Z  and 
G-Z  plots.   A  portion  of  the  actual  XBT  trace  is  shown  for  three  of  the 
stations.   In  these  three  cases,  there  is  excellent  visible  correlation 
between  the  actual  XBT  trace  and  the  reconstituted  T-Z  plots  which  shows 
that  MICROXBT  can  resolve  very  small  temperature  features.   In  no  case 
was  there  an  isothermal  layer  greater  than  140  feet  thick. 

2 .  Andaman  Sea 

In  the  Andaman  Sea,  near  the  Straits  of  Maalacca,  a  series  of 
15  XBT  stations  from  the  18  to  20  July  1969  period  were  used.   The 
stations  were  within  a  radius  of  12  miles.   Water  depths  varied  between 
290  and  380  feet.   This  oceanic  region  could  be  considered  somewhat 
similar  to  the  Atlantic  Coastal  region  above. 
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;ure  39.   T-Z  and  G-Z  Plots  for  Atlantic  Coastal  Region  6  and  7  December  1969 
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.re  40.   T-Z  and  G-Z  Plots  for  Atlantic  Coastal  Region  7 
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ure  41.   T-Z  and  6-Z  plots  for  Atlantic  Coastal  Region  9  and  10  December  1969 
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Figure  42.   T-Z  and  G-Z  Plots  for  Atlantic  Coastal  Region  11 


December  1969 
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The  T-Z  plots  (Figures  43  through  50)  contain  a  mixed  layer  to 
200  or  more  feet  except  for  one  observation  at  1200Z  on  the  19th  (Figure 
48b).   The  mixed  layer  transients  are  marked  with  spikes  in  the  G-Z  plots. 
Figure  50  shows  transients  with  the  greatest  magnitude  in  the  upper  100 
feet.   The  seasonal  thermocline  is  characterized  by  two  steps  (Figure 
44,  for  example).   The  G-Z  plots  in  most  cases  show  this  second  step, 
averaging  at  a  depth  about  300  feet,  with  good  resolution.   Some  strong 
gradients  are  found  in  the  200  to  300  feet  depth  region  in  every  G-Z 
plot  except  in  Figure  48b.   This  figure  appears  to  indicate  flow  over  a 
hump  or  have  erroneous  depth  values.   The  fall-rate  of  the  XBT  probe, 
appearing  erroneously  in  the  T-Z  plot,  indicates  a  possible  malfunction— 
of  the  probe  or  the  recorder  for  this  one  cast. 
3 .   Pacific  Oceanic  Region 

A  series  of  XBT  traces,  taken  over  a  week  period  in  October  1968, 
in  the  deep  waters  (over  100  nm  off  the  West  Coast  near  San  Diego, 
California)  were  used  for  this  observation.   A  composite  T-Z  plot  of 
the  actual  XBT  temperature-depth  traces  is  shown  in  Figure  51.   This 
region  is  sub-tropical  oceanic  [Tully,  1964].   Plots  of  T-Z  and  G-Z 
(Figures  52  through  65)  are  presented  from  the  surface  to  900  feet.   The- — 
seasonal  thermocline  is  evident  in  all  plots.   The  strong  thermocline  is 
clearly  shown.   Some  transients  are  notable.   Few  inversion  layers  were 
present,  and,  generally,  they  were  found  only  in  the  near-surface  layer. 
The  G-Z  plot  is  characterized  by  strong  gradients  in  the  seasonal  thermo- 
cline and  by  weak  gradients  in  the  main  thermocline,  each  separate  layer 
separated  by  isothermal  waters  from  one  to  forty  feet  thick. 

The  mixed  layer  contained  numerous  temperature  fluctuations  in 
most  cases.   Some  inversions,  in  particular  the  one  at  0600Z  on  14  October 
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Figure  43.   T-Z  and  G-Z  Plots  for  Andaman  Sea  Region  for  18  July  1969. 
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ure  44.   T-Z  and  G-Z  Plots  for  Andaman  Sea  Region  for  18  July  1969. 
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Figure  45.   T-Z  and  0-1  Plots  for  Andaman  Sea  Region  for  18  and  19  July  1969. 
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dre   46.   T-Z  and  G-Z  Plots  for  Andaman  Sea  Region  for  19  July  1969. 
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Figure  47.   T-Z  and  G-Z  plots  for  Andaman  Sea  Region  for  19  July  1969 
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.gure  48.   T-Z  and  G-Z  Plots  for  Andaman  Sea  Region  for  19  July  1969. 
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Figure  49.   T-A  and  G-Z  Plots  for  Andaman  Sea  Region  for  19  July  1969 
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Figure  50.   T-Z  and  G-Z  Plots  for  Andaman  Sea  Region  for  20  July  1969 
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Figure  51.   Composite  T-Z  Plot  from  Actual  XBT  Traces  Down  to  900  feet 
for  Region  100  run  off  San  Diego.   Data  are  given  in  time 
sequence  from  10  through  16  October  1968. 
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gure  52.   T-Z  and  GtZ  Plots  for  Eastern  North  Pacific  for  1800Z  10  October  1968, 
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Plgure  53.   T-Z  and  G-Z  PlotB  for  Eastern  North  Pacific  for  O0OOZ  11  October  1 
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gure  54.   T-Z  and  G-Z  Plots  for  Eastern  North  Pacific  for  0000Z  12  October  1968. 
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Figure  55.   T-Z  and  G-Z  Plots  for  Eastern  North  Pacific  for  0600Z  12  October  1968. 
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Figure  56.   T-Z  and  G-Z  Plots  for  Eastern  North  Pacific  for  1800Z  12  October  1968. 
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Figure  57.   T-Z  and  G-Z  Plots  for  Eastern  North  Pacific  for  0000Z  13  October  1968 
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Figure  58.   T-Z  and  G-Z  Plots  for  Eastern  North  Pacific  for  0600Z  13  October  1968, 
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Plgure  59i   T-Z  and  G-Z  Plots  for  Eastern  North  Pacific  for  1200Z  13  October  1968 
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igure  60.   T-Z  and  G-Z  Plots  for  Eastern  North  Pacific  for  1800Z  13  October  1968. 
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Figure  61.   T-Z  and  G-Z  Plots  for  Eastern  North  Pacific  for  0600Z  14  October 
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Lgure  62.   T-Z  and  G-Z  Plots  for  Eastern  North  Pacific  for  1200Z  14  October  1968, 
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Figure  63.   T-Z  and  G-Z  Plots  for  Eastern  North  Pacific  for  1800Z  14  October  1968 
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igure  64.   T-Z  and  G-Z  Plots  for  Eastern  North  Pacific  for  0000Z  15  October  1968 
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Figure  65.   T-Z  and  G-Z  Plots  for  Eastern  North  Pacific  for  0600Z  16  October  1968 
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1968  at  120  feet  was  clear  on  both  the  bathythermogram  and  the  T-Z  plot. 
The  steps  are  quite  irregular  on  the  T-Z  plot  but  not  as  irregular  as  on 
the  G-Z  plot. 

4.   Gulf  Stream  Region 

Thirteen  stations,  during  the  period  11-13  August  1969,  were 
analyzed  for  the  Gulf  Stream  region  (Figures  66  through  80).   The  ship's 
track  crossed  the  Gulf  Stream  boundary  several  times. 

a.  Interior  of  Gulf  Stream 

For  the  interior  of  the  Gulf  Stream  (Figures  70  and  71), 
the  mixed  layer  is  distinct  and  a  typical  seasonal  and  main  thermocline 
are  shown.   A  "boot"  in  the  100  foot  layer  is  clearly  shown  by  positive 
gradient  values  in  several  of  the  figures  and  the  actual  XBT  traces. 

b.  Frontal  Region 

Strong  gradients  are  the  rule  for  this  region  and  they  are 
quite  distinct  in  both  the  T-Z  and  G-Z  plots.   The  contrast  between 
the  frontal  and  interior  regions  of  the  Gulf  Stream  are  quite  apparent. 
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Figure   66.      T-Z  and  G-Z   Profile    for   Transects   across   Northern  Boundary   of   the   Gulf  Stream.      Actual  XBT  T-Z   Trace    Inset. 
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Figure  67.   T-Z  and  G-Z  Profile  for  Transects  across  Northern  Boundary  of  the  Gulf  Stream.   Actual  XBT 
Trace  Inset. 

103 


C2S    <xl0>F      00, 


Z(xlO   )£t 

Plgure   68.      T-Z  and  G-Z    Profile    for   Transects   across   Northern  Boundary   of   the  Gulf  Stream.      Actual  XBT  T-Z 
Trace    Inset. 
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Figure  69.   T-Z  and  G-Z  Profile  for  Transects  across  Northern  Boundary  of  the  Gulf  Stream.   Actual  XBT  T-Z 
Trace  Inset. 
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Figure  70.   T-Z  and  G-Z  Profile  for  Transects  across  Northern  Boundary  of  the  Gulf  Stream.   Actual  XBT  T-Z 
Trace  Inset. 
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Figure  71.   T-Z  and  G-Z  Profile  for  Transects  across  Northern  Boundary  of  the  Gulf  Stream.   Actual  XBT  T-Z 


Trace  Inset. 


107 


,  :  I  I  i  I  :  I  iTTI  i  J_r— 
' L.  J L-LJ-LI—J ;  —  6ri'\ 

_> i  /  T  '  i  i : ! ;;j.j..: 


c 


_  1 


120869 
3821N 


193; 

05244h 


Z(xl0  )ft 
Figure  72.   T-Z  and  "G^Z  Profile  for  Transects  across  Northern  Boundary  of  the  Gulf  Stream.   Actual  XBT  T-Z  Trace  Inset, 
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e  73.   T-Z  and  G-Z  Profile  for  Transecti~across  Northern  Boundary  of  the  Gulf  Stream.   Actual  XBT  T-Z  Trace  Inset. 
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Figure    74.      T-Z  and  G-Z   Profile    for   Tra 
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nsects  across  Northern.  Boundary  of  the  Gulf  Stream. 
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Figure  75.   T-Z  and  G-Z  Profile  for  Transects  across  Northern  Boundary  of  the  Gulf  Stream. 
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Figure    76.      T-Z  and   G-Z    Profile    for  Transects   across   Northern  Boundary   of   the   Gulf  Stream. 
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Figure  77.   T-Z  and  G-Z  Profile  for  Transects  across  Northern  Boundary  of  the  Gulf  Stream. 

113 


Z(xlO    )ft 
Figure   78-    f-Z  and  G-Z   Profile    for   Transects   across   Northern  Boundary   of   the   Gulf  Stream. 
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Figure   79.      T-Z  and   G-Z   Profile    for   Transects   across   Northern  Boundary   of   the   Gulf  Stream. 
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Figure  80.   T-Z  and  G-Z  Profile  for  Transects  across  Northern  Boundary  of  the  Gulf  Stream. 
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VI.   CONCLUSIONS 

The  results  give  a  measure  of  the  microthermal  structure  obtained 
with  an  XBT  by  the  NDM. 

A.   PHYSICAL  SIGNIFICANCE 

Microthermal  structure  statistics  have  not  been  adequately  explored. 
This  thesis  has  proposed  an  easy  method  to  observe  the  fine-scale  thermal 
structure.   Further,  this  proposed  objective  analysis  scheme  does  not  use 
the  T-Z  values  but  rather  the  G-Z  values.   It  is  the  G-Z  distribution 
which  is  important.   The  T-Z  can  always  be  reconstructed  if  one  additional 
parameter,  the  SST,  is  known. 

1.  Internal  Waves  -  Va'isa'la  Frequency 

The  Brunt-Vaisala  frequency  is  the  maximum  frequency  of  oscil- 
lation of  free  internal  waves.   It  may  be  calculated  from  the  density 
gradient,  or  from  the  parameters  which  define  density  gradient.   In 
isohaline  water,  the  NDM  gives  a  direct  measure  of  density  (or  stability). 
In  such  waters  it  could  be  used  to  set  an  upper  limit  on  the  maximum 
frequency  of  oscillation  of  internal  waves.   Further,  a  time  series 
study  would  clearly  illustrate  the  amplitudes  likely  to  be  encountered 
since  the  sheets  indicated  on  the  G-Z  diagram  are  readily  identified. 
A  correlation  study  of  G-Z,  Brunt-Vaisala  versus  depth  diagrams  would 
be  useful. 

2 .  Oceanic  Fronts 

This  difficult  region  to  analyze  microscopically  could  possible 
be  better  explored  by  using  the  NDM  for  the  microthermal  processes.   The 
NDM  explores  the  continuously  changing  conditions  in  frontal  regions. 
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Frontal  regions  are  usually  indicated  by  extremely  complex  G-Z  plots 
which  illustrate  the  extreme  gradients,  both  positive  and  negative, 
encountered.   Based  on  the  admittedly  small  number  of  cases  examined, 
there  appears  to  be  a  tendency  of  the  layers  to  maintain  their  character- 
istics and  slide  by,  over  or  under  each  other  rather  than  mix. 
3 .   Sound  Propagation 

The  effect  of  microthermal  structure  on  sound  propagation 
require  better  statistical  knowledge  of  the  oceans.   The  NDM  gives 
emphasis  to  the  strongest  gradient  regions. 

The  microthermal  effects  on  index  or  refraction  variations 
(density  fluctuations)  must  be  determined.   Using  the  NDM  with  salinity 
variations  (from  an  STD  cast)  and  the  NDM  with  sound  velocity  variations 
(from  an  SVTP  cast)  taken  from  the  two  instruments  simultaneously,  cross- 
correlation  techniques  could  shed  some  light  on  the  unresolved  micro- 
processes  in  the  ocean.   In  particular,  what  are  the  exact  physical  and 
theoretical  correlations  among  the  salinity,  temperature,  and  sound 
velocity . 

B.   SUGGESTED  STUDIES 

Possible  studies  for  the  future  have  suggested  themselves  while 
developing  this  thesis. 

1 .   Mathematical  Models 

Various  mathematical  models  attempted  could  be  explored  further 
because  each  method  had  some  merit.   Suggested  models  are  limited  by 
their  utility. 

a.   Analysis  Methods 

The  microthermal  structure  recorded  by  XBTs  can  be  looked  at 
as  a  "roughness"  factor,  the  noise  level  of  the  water  structure.   Further 
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work  is  required  to  determine  such  a  factor  for  practical  operational  use. 
This  roughness  factor  should  be  computed  as  the  XBT  trace  is  being 
analyzed  at  FNWC  by  a  program  like  DIGITBT.   This  factor  could  be  one 
number  to  describe  objectively  the  microthermal  structure  variations  on 
the  G-Z  trace . 

Variations  in  the  NDM  should  be  used  with  the  variables  that 
have  been  arbitrarily  set.   For  example,  the  depth  interval  of  AT/ft 
at  one  foot  should  be  varied  to  find  the  effects  a  larger  or  smaller 
depth  interval  would  have  on  the  G-Z  profiles.   As  the  steps  become 
smaller,  the  exact  derivative  is  approached.   Larger  steps  should  smooth 
the  spikes  shown  in  the  G-Z  figures.   In  some  cases,  this  smoothed  G-Z 
profile  might  be  described  with  curve  fitting  techniques  as  suggested  by 
Boston  [1966]  and  Grosfils  [1968]. 

MICROXBT's  digitizing  interval  of  0.01  inches  should  be 
changed  to  0.005  inches  and  the  program  modified  to  use  this  interval. 
This  would  yield  more  detailed  thermal  structure  than  MICROXBT.   Obviously, 
one  cannot  make  the  interval  finer  than  the  noise  level  of  the  sensor  or 
recording  instrument. 

b.   Physical  Processes 

The  G-Z  profiles  from  the  NDM  emphasize  the  changing  physical 
features  which  are  obscure  in  the  T-Z  profiles.   By  using  G-Z  profiles 
for  parameters  other  than  temperature,  clarity  in  determining  distri- 
butions of  chemical  parameters  would  be  similarly  emphasized.   Limitations 
in  use  of  the  NDM  on  other  chemical  properties  should  be  further 
investigated. 

Water  mass  identification  is  possible  with  G-Z  profiles. 
Additional  work  is  required  in  this  area. 
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2  .   Experimental  Work 

A  microthermal  sheet  probe  could  be  designed  for  investigating  the 
chemical  and  physical  processes  inside  the  sheet  and  inside  the  layers. 
A  comparison  of  actual  physical  property  differences  should  be  made. 

A  modified  expendable  temperature  gradient  (XTG)  meter  could  be 
designed  using  two  thermistors  to  obtain  the  temperature  differences  in 
the  probe.   The  standard  XBT  system  equipment  already  installed  could  be 
changed  to  give  analog  G-Z  profiles  to  the  operational  user  without  any 
de lay . 

a.  Time  Series 

The  XTG  meter  could  be  used  to  construct  G-Z  variability  at 
one  location  by  conducting  a  time  series  analysis.   This  needs  to  be 
done  in  order  to  obtain  a  roughness  factor,  an  objective  microstructure 
measuring  parameter. 

b.  Transect 

An  XTG  meter  could  be  used  in  transects  across  various  water 
masses  to  obtain  the  definite  characteristics  for  the  water  masses.   The 
XTG  could  be  used  to  follow  interesting  features  on  the  G-Z  plots.   More 
important,  however,  is  that  the  spacial  characteristics  of  microthermal 
structure  could  easily  be  determined. 
3 .   Handling  of  Data 

When  the  microthermal  structure  roughness  factor  is  determined, 
this  factor  should  be  made  an  integral  part  of  the  records  for  that  XBT 
trace.   This  information  should  be  included  as  part  of  the  ID  data. 

Aboard  Navy  ships  the  AN/UNQ-7  tape  recorder  could  be  used  to 
obtain  magnetic  tape  records  of  the  analog  XBT  or  XTG  signals.  These 
tapes  could  be  submitted  with  the  traces  to  FNWC  or  NPS  for  further  analys 
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4 .   Naval  Applications 

Results  from  microthermal  studies  show  that  high  frequency  sonar 
performance  as  influenced  by  the  microstructure  receive  the  greatest 
microstructure  effects.   More  knowledge  of  the  microstructure  is  needed 
before  the  torpedo  and  mine-hunting  sonar  performance  could  be  improved. 
With  a  measure  of  the  microstructure,  the  performance  of  a  sonar  system 
could  be  predetermined. 

Tens  of  thousands  of  XBT  traces  are  on  file  at  FNWC.   New  data 
are  arriving  daily.   This  thesis  is  an  attempt  to  show  how  this  data  can 
be  used  both  operationally  and  scientifically.   Obviously  it  represents 
only  a  beginning  in  the  way  of  microstructure  analysis  and  operational 
utilization. 
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APPENDIX  A 


MICROXBT  FORTRAN   IV   PROGRAM  FOR  CDC    6500 
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//(Job  name) 
//  EXEC 
//F0RT.SYSIN 


/* 

//G0.FTO6FOO1 

//G0.FTO4FOO1 

II 

II 

//G0.FTO4FOO2 

II 

II 

//G0.SYSIN 


J0B  (Acct.  no.,  Job  no.,  Section),  '(Your  name)'  Note  1 
F0RTCLGP,REGION=1OOK,TIME.G0=2  Note  2 

DD  * 

FORTRAN  Source  Deck 

Note  3 
DD  SYS0UT=A,SPACE(CYL,6)  Note  4 

DD  UNIT=24OO-1,V0L=SER=NPSXXX,IABEL=(,NL),  Note  5 
DCB=(DEN=1,RECFM=F,LREC1>136,BLKSIZE=136,TRTCH=ET),  "  5 
DISP=(0LD,KEEP)  Note  5 

DD  AFT=FTO4FOO1,V0L=SER=NPSXXX,IABEL(2,NL),  Note  6 
DCB=(DEN=1,RECFM=F,LRECL=136,BLKSIZE=136,TRTCH=ET),  "  6 
DISP=(0LD,KEEP)  Note  6 


DD  * 


Note  7 


DATA  Cards,  if  any 
/*  (Orange  Card) 
Notes 

1.  Standard  J0B,  green  card,  for  program  submission. 

2.  If  not  plotting  or  card  punching,  the  F0RTCLGP  can  be  changed 
to  F0RTCLG. 

3.  Must  be  a  white  card. 

4.  Optional  card  for  additional  storage  for  output. 

5.  The  XXX  following  NPS  is  the  NPS  tape  number.   These  three  cards 
are  for  reading  the  first  file  on  the  input  tape.   These  cards 
assign  input  unit  4  for  the  source  programs  READ  statements. 
Further  note  that  the  record  length  is  136  characters,  all  of 
which  may  be  read.   This  is  the  standard  record  length  from  the 
CDC  3600. 

6.  For  a  third  file  to  be  read,  these  three  cards  must  be  repeated 
with  F002  changed  to  F003  and  IABEL(2,NL)  changed  to  1ABEL(3,NL) 
and  so  forth  for  additional  files.   Also  the  UNIT=2400-1  should 
be  changed  to  AFT=FT04F001  for  all  files  after  the  first  file. 
If  only  one  file  is  to  be  read,  then  these  three  cards  are  not 
required. 

7.  Not  required  if  there  are  no  input  data  cards. 
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//(Job  name)    J0B  (Acct.  no.,  Job  no.,  Section),  '(Your  name)'      Note  1 

//DUMPTAPE     EXEC  PGM=IEBPTPCH 

//SYS PRINT     DD  SYS0UT=A 

//SYSUT1       DD  UNIT=24OO-1,DISB=0LD,IABE]>(,NL),V0LUME=SER=NPSXX,  Note  2 

//  DCB=(DEN=1,RECFM=F,LRECD=136,BLKSIZE=136,TRTCH=ET)    Note  2 

//SYSUT2       DD  SYS0UT=A 

//SYSIN        DD  *,DCB=BLKSIZE=80 

PRINT  MAXFLDS=1,ST0PAFT=1OOO  Note  3 

RECORD  FIELDS=(120)  Note  4 

/*  (Orange  Card) 

Notes 

1.  Standard  J0B,  green  card,  for  program  submission. 

2.  If  more  than  one  file  is  on  the  tape  then  separate  two  cards  for 
each  file  identified  by  file  number  by  the  file  number  in  the  IABEL 
(X,NL)  where  S  is  the  file  number,  which  is  implied  1  when  omitted. 

3.  The  ST0PAFT  limits  the  number  of  lines  of  printed  output,  in  this 
example  1000  lines  is  used  before  printing  is  terminated. 

4.  The  field  is  recorded  at  the  maximum  amount  that  this  program  will 
dump  which  is  120  characters.   If  there  are  any  characters  between 
120  and  133  in  the  records  then  this  routine  will  not  dump  those 
characters . 


145 


BIBLIOGRAPHY 


Arthur  D.  Little  Inc.,  1965.   Experimental  Evaluation  of  Expendable 
Bathythermographs ,  Report  4071165. 

Arthur  D.  Little  Inc.,  1966.   Internal  Waves :   Their  Influence  Upon 
Nava 1  Operat ions ,  Report  4090266. 

Arthur  D.  Little  Inc.,  1966a.   Expendable  Bathythermograph  (XBT)  System 
Evaluation  for  Tactical  Sonar  Application,  Report  4150866. 

Arthur  D.  Little  Inc.,  1966b.   Acoustic  Transmission  in  an  Ocean  Surface 
Duct,  Report  4171166. 

Arthur  D.  Little  Inc.,  1967.   Internal  Waves  in  the  Ocean,  Report  4321267 

Arthur  D.  Little  Inc.,  1968.   An  Assessment  of  Expendable  Bathythermo- 
graph (XBT)  Data  and  Recorder  Requirements,  Report  4250268. 

Arthur  D.  Little  Inc.,  1968a.   Experimental  Evaluation  of  Deep  (6000  ft.) 
Expendable  Bathythermographs,  Report  4391168. 

Ayers ,  E.,  Wolff,  P.  M.  ,  Carstensen,  L.  ,  and  Ayers ,  H.  C,  1966.   A  Ray 
Tracing  Program  for  a  Digital  Computer,  paper  presented  at  NATO 
Symposium  on  the  Operational  Evaluation  of  ASW  Weapons  Systems, 
Paris,  France  on  28  June  1966. 

Barakos ,  P.  A.,  1968.   Acoustic  Wave  Scattering  and  Refraction  by 

Internal  Waves  in  a  Model  Tank,  U.  S.  Navy  Underwater  Sound  Labora- 
tory Report  905,  Ph.D.  Thesis,  University  of  Washington,  Seattle. 

Bauer,  R.  A.,  1969.   Fleet  Numerical  Weather  Central's  XBT  Digitizing 
System,  Computer  Applications  Report. 

Berkerle,  J.  C,  Reitzell,  J.  S.,  Payne,  R.  E.,  and  Hersey,  J.  B.,  1968. 

Variations  in  Sound  Between  Bermuda  and  the  Antilles,  lectures  given 

at  AGU  meeting  on  20  April  1966,  Woods  Hole  Oceanographic  Institute 
Reference  68-19  (Unpublished  manuscript). 

Bergmann,  P.  G.,  1946.   "Propagation  of  Radiation  in  a  Medium  with 
Random  Inhomogeneities ,"  Physical  Review,  v.  70,  p.  486-492. 

Boston,  N.  E.  J.,  1966.   Objective  Definition  of  the  Thermocline,  Texas 
A&M  Reference  66-19T,  Project  286-D. 

Choate ,  T.  V.,  1970.   (Marketing  Manager,  Sippican  Corporation),  personal 
communication . 

Cohen,  J.  S.,  and  Weinberg,  H.,  1969.   CONGRATS  Ray  Plot  Program,  USL 
Technical  Memorandum  2070-110-69. 


146 


Cook,  J.  C,  and  Kenyon,  K.  E.,  1963.   "Fast-Response  Thermistor  Probes 
for  Temperature  Microstructure  Studies  at  Sea,"  Review  of  Scientific 
Instruments ,  v.  34,  p.  496-499. 

Cooper,  J.  W. ,  and  Stommel,  H. ,  1968.   "Regularly  Spaced  Steps  in  the 
Main  Thermocline  near  Bermuda,"  Journal  of  Geophysical  Research, 
v.  73,  p.  5849-4854. 

Cooper,  L.  H.  N. ,  1967.   "Stratification  in  the  Deep  Ocean,"  Sci.  Prog. 
(Oxf .),  v.  55,  p.  73-90. 

Demeo,  R.  P.,  1968.   The  Validity  of  Expendable  Bathythermograph  Mea- 
surements ,  Paper  presented  at  MTS  Meeting,  4th,  Miami,  Florida. 

Denner,  W.  W.  ,  1966.   A  Comparison  of  the  Sippican  and  General  Motors 
Expendable  Bathythermograph  to  a  Mechanical  Bathythermograph,  Fleet 
Numerical  Weather  Central  Technical  Note  16. 

Denner,  W.  W. ,  1968.   "Separation  of  the  Residual  Instrument  Noise  from 
the  Significant  Variability  for  Expendable  Bathythermographs," 
Marine  Sciences  Instrumentation,  v.  4,  Plenum  Press,  p.  635-641. 

Denner,  W.  W.,  1969.   The  T-S  Gradient  Method,  a  New  Method  of  Computing 
Geostrophic  Currents  over  Large  Ocean  Areas,  U.  S.  Naval  Postgraduate 
School  Technical  Report  NPS-58DW9072A . 

General  Motors  Defense  Laboratories,  1965.   Observations  on  Expendable 
Instrumentation. 

DiNapoli,  F.  R. ,  1969.   Acoustic  Propagation  in  a  Stratified  Medium, 
USL  Report  1046,  Ph.D.  Thesis,  University  of  Rhode  Island. 

Einstein,  L.  T. ,  and  Cohen,  J.  S.,  1969.   Acoustic  Field  Modeling  with 
Constant  and  Continuous  Velocity  Gradients,  USL  Technical  Memorandum 
2070-237-69. 

Frankignoul,  C,  1969.   Note  on  Internal  Waves  in  a  Simple  Thermocline 
Model,  Woods  Hole  Ref.  69-47  (Unpublished  manuscript). 

Glennon,  A.  N. ,  1965.   "Practicalities  of  Sonar  Range  Prediction,"  U.  S. 
Naval  Institute  Proceedings,  v.  91,  p.  128-132,  March  1965. 

Gouzie,  M.  W. ,  Sanders,  M.  R. ,  and  Littlehale,  A.  D.,  1966.   Laboratory 
and  Engineering  Evaluation  of  Three  Models  of  Expendable  Bathythermo- 
graphy,  USL  Report  732. 

Grant,  H.  L.  ,  Moilliet,  A.,  and  Vogel,  W.  M. ,  1968.   "Some  Observations 
of  the  Occurrence  of  Turbulence  in  and  above  the  Thermocline," 
Journal  of  Fluid  Mechanics,  v.  34,  p.  443-448. 

Green,  E.  L.  ,  1967.   Diffraction  Theory  of  Lens  less  Correlation,  USL 
Report  842. 


147 


Grosfils,  E.  F.,  1968.   Objective  Digital  Analysis  of  Bathythermograph 

Traces ,  M.S.  Thesis,  Naval  Postgraduate  School,  Monterey,  California. 

Haurwitz,  B.,  Stommel,  H. ,  and  Munk,  W.  H. ,  1959.   "On  the  Thermal  Unrest 
in  the  Ocean,"  in  The  Atmosphere  and  the  Sea  in  Motion,  ed.  by  B . 
Boiln,  Rockefeller  Institute  Press,  N.  Y.,  p.  74-94. 

Hazelworth,  J.  B.,  1966.   Quantitative  Analysis  of  Some  Bathythermograph 
Errors,  U.  S.  Naval  Oceanographic  Office  TR-180,  ASWEPS  Report  111. 

Hecht,  A.,  and  White,  R.  A.,  1968.   "Temperature  Fluctuations  in  the 
Upper  Layer  of  the  Ocean,"  Deep-Sea  Research,  v.  15,  p.  339-353. 

Hubert,  W.  E.,  and  T.  Laevastu.,  1966.   Short -Period  Changes  and 

Anomalies  of  Temperature  in  the  Oceans  and  Their  Effects  on  Sound 
Propagation,  Fleet  Numerical  Weather  Facility  Technical  Note  10. 

Ingham,  M.  C,  1968.   "The  'Mixed  Layer'  in  the  Western  Tropical  Atlantic 
Ocean,"  Bulletin  of  Marine  Science,  v.  18,  p.  561-571. 

James,  R.  W. ,  1966.   Ocean  Thermal  Structure  Forecasting,  NavOceanO, 
SP-105,  ASWEPS  Manual  v.  5. 

James,  R.  W. ,  1967.   Data  Requirements  for  Synoptic  Sea  Surface  Tempera- 
ture Analyses  ,  NavOceanO,  SP-111. 

Knollman,  G.  C,  1964.   "Wave  Propagation  in  a  Medium  with  Random, 

Spheroidal  Inhomogeneities ,"  Journal  of  the  Acoustical  Society  of 
America,  v.  36,  p.  681-688. 

Krauss ,  W. ,  1966.  Internal  Tides  off  Southern  California,  U.  S.  Naval 
Electronics  Laboratory  Report  1389. 

LaFond,  E.  C,  1959.  Slicks  and  Temperature  Structure  in  the  Sea,  NEL 
Report  937. 

LaFond,  E.  C,  and  Cox,  S.  S.,  1962.   "Internal  Waves,"  in  The  Sea,  v.  1, 
ed.  by  W.  M.  Hill,  Interscience  Publishers,  p.  731-763. 

LaFond,  E.  C,  and  LaFond,  K.  G.  ,  1966.   Vertical  and  Horizontal 
Structures  in  the  Sea,  NEL  Report  1395. 

Laevastu,  T.,  1960.  "Factors  Affecting  the  Temperature  of  the  Surface 
Layer  of  the  Sea,"  So'c.  Scient.  Fennica.  Comm.  Phys-Mathem,  v.  25, 
p.  1-136. 

Lee,  0.  S.,  1961.  "Effect  of  an  Internal  Wave  on  Sound  in  the  Ocean," 
Journal  of  the  Acoustical  Society  of  America,  v.  33,  p.  677-681. 

Lee,  0.  S.,  1961a.  "Observations  on  Internal  Waves  in  Shallow  Water," 
Limnology  and  Oceanography,  v.  6,  p.  312-321. 


148 


Liebermann,  L.  N. ,  1951.   "The  Effect  of  Temperature  Inhomogeneities 

in  the  Ocean  on  the  Propagation  of  Sound,"  Journal  of  the.  Acoustical 
Society  of  America,  v.  23,  p.  563-57-. 

Lovett,  J.  R. ,  1968.   "Vertical  Temperature  Gradient  Variations  Related 
to  Current  Shear  and  Turbulence,"  Limnology  and  Oceanography,  v.  13, 
p.  127-142. 

Meyer,  R.  F. ,  and  Romberg,  B.  W. ,  1963.   Acoustic  Scattering  in  the  Ocean, 
Arthur  D.  Little  Report  1360863. 

Miles,  J.  W.  ,  1961.   "On  the  Stability  of  Heterogeneous  Shear  Flows," 
Journal  of  Fluid  Mechanics,  v.  10,  p.  496-501. 

Miles,  J.  W.  ,  1963.   "On  the  Stability  of  Heterogeneous  Shear  Flows, 
Part  2,"  Journal  of  Fluid  Mechanics,  v.  16,  p.  209. 

Mintzer,  D.,  1953.   "Wave  Propagation  in  a  Randomly  Inhomogeneous  Medium, 
I,"  Journal  of  the  Acoustical  Society  of  America,  v.  25,  p.  922-927. 

Mintzer,  D.,  1953a.   "Wave  Propagation  in  a  Randomly  Inhomogeneous 

Medium,  II,"  Journal  of  the  Acoustical  Society  of  America,  v.  25, 
p.  1007-1111. 

Mintzer,  D.,  1954.   "Wave  Propagation  in  a  Randomly  Inhomogeneous  Medium, 
III,"  Journal  of  the  Acoustical  Society  of  America,  v.  26,  p.  186-190. 

Mooers ,  C.  N.  K. ,  1966.   A  Study  of  Ocean  Variability  at  a  Depth  of  One 
Mile,  USL  Report  693. 

National  Oceanographic  Data  Center,  1968.   Unclassified  letter  to  all 
Sippican  XBT  Users,  Subject:   Procedures  for  Disposition  of  XBT 
Traces,  21  February  1968. 

National  Oceanographic  Data  Center,  1968a.   Instructions  for  Preparing 
the  Bathythermograph  Log  Sheet,  NODC  EXP-3167/10-A  (Rev.  1-68). 

National  Defense  Research  Committee,  1946.   The  Application  of  Ocean- 


ography to  Subsurface  Warfare,  Summary  Technical  Report  of  Division 
6 ,  v .  6A . 

Naval  Weather  Service,  1967.   U.  S.  Naval  Weather  Service  Computer 

Products  Manual,  Naval  Air  Systems  Command  Report  NAVAIR  50-1G-522. 

Neshyba,  S.,  Neal,  V.  T. ,  and  Denner,  W.  W. ,  1969.   "The  Significance 
of  Temperature  Stratification  in  the  Arctic,"  paper  presented  at 
Symposium  on  Military  Oceanography,  6th,  Seattle,  Washington. 

Officer,  C.  B.,  1958.   Introduction  to  the  Theory  of  Sound  Transmission 
with  Application  to  the  Ocean,  McGraw-Hill. 

Osborn,  T.  R. ,  1969.   Oceanic  Fine  Structure,  Ph.D.  Dissertation, 
University  of  California  at  San  Diego. 


149 


Pao,  Y-H. ,  1969.   "Origin  and  Structure  of  Turbulence  in  Stably  Strat- 
ified Media,"  Clear  Air  Turbulence  and  Its  Detection,  Plenum  Press, 
p.  73-99. 

Pederson,  M.  A.,  1961.   "Acoustic  Intensity  Anomalies  Introduced  by 
Constant  Velocity  Gradients,  Journal  of  the  Acoustical  Society  of 
America,  v.  33,  p.  465-474. 

Pedersen,  J.  A.,  and  Gordon,  D.  F.,  1965.   "Normal-Mode  Theory  Applied 
to  Short-Range  Propagation  in  an  Underwater  Acoustic  Surface  Duct," 
Journal  of  the  Acoustical  Society  of  America,  v.  37,  p.  105-110. 

Pipp,  A.  T.,  1964.   "Fine  Structure  and  Stability  of  the  Sound  Channel 
in  the  Ocean,"  Journal  of  the  Acoustical  Society  of  America,  v.  36, 
p.  1948-1953. 

Pingree ,  R.  D.,  1969.   "Small-Scale  Structure  of  Temperature  and  Salinity 
Near  Station  Cavall,"  Deep-Sea  Research,  v.  16,  p.  275-295. 

Pond,  S.,  1968.   "Microscale  Phenomena  in  the  Ocean  and  the  Atmospheric 
Boundary  Layer,"  Proceedings  of  Region  Six  IEEE  Conference  l-C-4, 
p.  1-10. 

Potter,  D.  S.,  and  Murphy,  S.  R. ,  1957.   "On  Wave  Propagation  in  a 

Random  Inhomogeneous  Medium,"  Journal  of  the  Acoustical  Society  of 
America,  v.  29,  p.  197-198. 

Rasmussen,  R.  A.,  1962.   "Application  of  Thermistors  to  Measurements  in 
Moving  Fluids,"  Review  of  Scientific  Instruments,  v.  33,  p.  38-42. 

Rattray,  M. ,  1962.   "Interpolation  Errors  and  Oceanographic  Sampling," 
Deep-Sea  Research,  v.  9,  p.  25-37. 

Reid,  W.  L. ,  1964.   Bathythermograph  Evaluation,  (Unpublished  manuscript). 

Richardson,  W.  M.,  1958.   Measurement  of  Thermal  Microstructure ,  Woods 
Hole  Ref  58-11  (Unpublished  manuscript). 

Roden,  G.  I.,  1963.   "On  Sea  Level,  Temperature,  and  Salinity  Variations 
in  the  Central  Tropical  Pacific  and  on  Pacific  Ocean  Islands," 
Journal  of  Geophysical  Research,  v.  68,  p.  455-472. 

Rossby,  C.  G.,  and  Montgomery,  R.  B. ,  1934.   "The  Layer  of  Frictional 

Influence  on  Wind  and  Ocean  Currents,"  Papers  in  Physical  Oceanography 
and  Meteorology,  v.  3. 

Sagar,  F.  H. ,  1960.  "Acoustic  Intensity  Fluctuations  and  Temperature 
Microstructure  in  the  Sea,"  Journal  of  the  Accoustical  Society  of 
America,  v.  32,  p.  112-121. 

Schulkin,  M. ,  1969.   The  Propagation  of  Sound  in  Imperfect  Ocean  Surface 
Ducts,  USL  Report  1013. 


150 


Sederowitz ,  W.  J.,  and  Favret ,  A.  G. ,  1969.  "Covariance  Between  Acoustic 
Signals  Traveling  Diverging  Paths  in  a  Random  Medium,"  Journal  of  the 
Acoustical  Society  of  America,  v.  45,  p.  386-391. 

Sheehy,  M.  J.,  1950.   "Transmission  of  24-Kc  Underwater  Sound  from  a 
Deep  Source,"  Journal  of  the  Acoustical  Society  of  America,  v.  22, 
p.  24-28. 

Shonting,  D.  H. ,  and  Kadis,  A.  L.  ,  1968.   The  Thermiprobe  :   A  System 

for  Measuring  Thermal  Microstructure  in  the  Sea,  paper  presented  at 
National  ISA,  Marine  Sciences  Instrumentation  Symposium,  4th, 
Pittsburg,  Pennsylvania  on  22  January  1968. 

Siedler,  G.,  1969.   "On  the  Fine  Structure  of  Density  and  Current  Dis- 
tribution and  its  Short-Time  Variations  if  Different  Areas,"  in 
Progress  in  Oceanography,  v.  5,  ed .  by  M.  Sears,  Pergmanon  Press, 
p.  81-94. 

Sippican  Corporation,  1968.   Instructions  for  Installation,  Operation 
and  Maintenance  of  Sippican  Expendable  Bathythermograph  System, 
Sippican  Report  R-467B. 

Skudrzyk,  E.,  1957.   "Scattering  in  an  Inhomogeneous  Medium,"  Journal  of 
the  Acoustical  Society  of  America,  v.  29,  p.  50-60. 

Smith,  E.  L. ,  1967.   A  Predictive  Horizontal-Temperature  Gradient  Model 
of  the  Upper  750  Feet  of  the  Ocean,  NEL  Report  1445. 

Spilhaus,  A.  F.,  1938.   "A  Bathythermograph,"  Journal  of  Geophysical 
Research,  v.  1,  p.  95-100. 

Stern,  M.  E.,  1968.   "T-S  Gradients  on  the  Micro-Scale,"  Deep-Sea 
Research,  v.  15,  p.  245-250. 

Stern,  M.  E.,  and  Turner,  J.  S.,  1969.   "Salt  Fingers  and  Convecting 
Layers,"  Deep-Sea  Research,  v.  16,  p.  495-511. 

Stommel,  H. ,  1966.   The  Gulf  Stream,  University  of  California  Press. 

Stommel,  H. ,  and  Federov,  K.  N.  ,  1967.   "Small  Scale  Structure  in 

Temperature  and  Salinity  Near  Timor  and  Mindanao,"  Tellus ,  v.  19, 
p.  306-325. 

Stone,  R.  G. ,  and  Mintzer,  D. ,  1962.   "Range  Dependence  of  Acoustic 
Fluctuations  in  a  Randomly  Inhomogeneous  Medium,"  Journal  of  the 
Acoustical  Society  of  America,  v.  34,  p.  647-653. 

Sverdrup,  H.  U.,  Johnson,  M.  W. ,  and  Fleming,  R.  H. ,  1942.   The  Oceans, 
Prentice-Hall. 

Tait,  R.  I.,  and  Howe,  M.  R. ,  1968.   "Some  Observations  on  Thermo-Haline 
Stratification  in  the  Deep  Ocean,"  Deep-Sea  Research,  v.  15,  p. 
275-280. 


151 


Tolstoy,  I.,  and  Clay,  C.  S.,  1966.   Ocean  Acoustics,  McGraw-Hill. 

Tully,  J.  P.,  1964.   "Oceanographic  Regions  and  Assessment  of  Tempera- 
ture Structure  in  the  Seasonal  Zone  of  the  Pacific  Ocean,"  J.  Fish. 
Res.  Bd.  Canada,  v.  21,  p.  941-970. 

Tully,  J.  P.,  1964a.  "Oceanographic  Regions  and  Processes  in  the  North 
Pacific  Ocean,"  in  Progress  in  Oceanography,  v.  2,  ed.  by  M.  Sears, 
Pergmanon  Press,  p.  68-84. 

Turner,  J.  S.,  1967.   "Salt  Fingers  Across  a  Density  Interface,"  Deep-Sea 
Research,  v.  14,  p.  599-611. 

Turner,  J.  S.,  and  Stommel,  H.,  1964.   "A  New  Case  of  Convection  in  the 
Presence  of  Combined  Vertical  Salinity  and  Temperature  Gradients," 
Proc .  Nat.  Academy  of  Science,  v.  52,  p.  49-53. 

Urick,  R.  J.,  1967.   Principles  of  Underwater  Sound  for  Engineers, 
McGraw-Hill,  p.  122-160. 

Urick,  R.  J.,  1967a.   Observations  of  the  Sound  Field  in  the  Deep  Sea, 
Naval  Ordnance  Laboratory  Technical  Report  67-103. 

Urick,  R.  J.,  and  Searf oss ,  C.  W. ,  1948.   The  Microthermal  Structure 
of  the  Ocean  Near  Key  West,  Florida,  Part  I,  Discussion,  Naval 
Research  Laboratory  Report  S-3392. 

Urick,  R.  J.,  and  Searf oss ,  C.  W. ,  1949.   The  Microthermal  Structure 
of  the  Ocean  Near  Key  West,  Florida,  Part  II,  Analysis,  Naval 
Research  Laboratory  Report  S-3444. 

Van  Vliet,  C.  J.,  1969.   An  Algorithm  for  Estimating  the  Isothermal 
Layer  Depth  of  a  Digitized  Bathythermograph,  NUC  TN  269. 

Weinberg,  H. ,  1967.   A  Curve  Fitting  Technique  for  Acoustic  Ray 
Analysis ,  USL  Report  805. 

White,  R.  A.,  1967.   "The  Vertical  Structure  of  Temperature  Fluctuations 
within  an  Oceanic  Thermocline ,"  Deep-Sea  Research,  v.  14,  p.  613-623. 

Whitemarsh,  D.  C,  Skrudzyk,  D.,  and  Urick,  R.  J.,  1957.   "Forward 
Scattering  of  Sound  in  the  Sea  and  Its  Correlation  with  the 
Temperature  Microstructure ,"  Journal  of  the  Acoustical  Society  of 
America,  v.  29,  p.  1124-1143. 

Wilson,  D.  F.,  Ricalzone,  L.  C,  and  Beckett,  R.  C,  1968.   Some 

Investigations  of  Sound-Scattering  Layers  Near  Key  West,  Florida, 
Naval  Research  Laboratory  Report  6682. 

Wolff,  P.  M.,  Tatro,  P.  R. ,  and  Megehee,  L.  D.,  1967(7).   The  FNWF 
Sound  Map  Program  (Unpublished  manuscript). 


152 


Woods,  J.  D. ,  1968.   "Wave-Induced  Shear  Instability  in  the  Summer 
Thermocline,"  Journal  of  Fluid  Mechanics,  v.  32,  p.  791-800. 

Woods,  J.  D.,  1968a.   "An  Investigation  of  Some  Physical  Processes 

Associated  with  the  Vertical  Flow  of  Heat  Through  the  Upper  Ocean," 
Meteorological  Magazine,  v.  98,  p.  65-73. 


153 


INITIAL  DISTRIBUTION   LIST 


No.    Copies 


1.  The  Oceanographer  of  the  Navy  1 
The  Madison  Building 

732  N.  Washington  Street 
Alexandria,  Virginia  22314 

2.  Commanding  Officer  1 
Fleet  Numerical  Weather  Central 

Monterey,  California  93940 

3.  Defense  Documentation  Center  20 
Cameron  Station 

Alexandria,  Virginia  22314 

4.  Library,  Code  0212  2 
Naval  Postgraduate  School 

Monterey,  California  93940 

5.  Professor  W.  W.  Denner,  Code  58Dw  4 
Department  of  Oceanography 

Naval  Postgraduate  School 
Monterey,  California  93940 

r  [ 

6.  Professor  N.  E.  J.  Boston,  Code  58Bb  1 
Department  of  Oceanography 

Naval  Postgraduate  School 
Monterey,  California  93940 

7.  Professor  T.  D.  Burnett,  Code  55Bu  1 
Department  of  Operations  Analysis 

Naval  Postgraduate  School 
Monterey,  California  93940 

8.  Department  of  Oceanography  3 
Naval  Postgraduate  School 

Monterey,  California  93940 

9.  LCDR  P.  M.  Magruder,  USN  4 
C OMC RUDE S PLOT  TEN 

Fleet  Post  Office 

New  York,  New  York  09501 


154 


Security  Classification 


DOCUMENT  CONTROL  DATA  -R&D 

(Security  classification  ol  title,  body  of  abstract  and  indexing  annotation  must  be  entered  when  the  overall  report  Is  classilied) 


irigin  a  ting   activity  (Corporate  author) 

Naval   Postgraduate   School 
Monterey,   California   93940 


2a.   REPORT    SECURITY    CLASSIFICATION 


Unclassified 


2b.    GROUP 


1EPORT    Tl  TLE 


Some  Characteristics  of  Temperature  Microstructure  in  the  Ocean 


descriptive  NOTES  (Type  ot  report  and. inclusive  dates) 

Master's   Thesis;  April    1970 

ku  thorisi  (First  name,  middle  initial,  last  name) 


Petyon  Marshall  Magruder,  Jr 


IEPORT    DATE 


April   1970 


7a.    TOTAL    NO.    OF    PAGES 


155 


7b.    NO.    OF    REFS 


109 


CONTRACT    OR    GRANT    NO. 


PROJEC  T    NO. 


•a.    ORIGINATOR'S    REPORT    NUMBER(S) 


9b.   OTHER  REPORT  NO(S)  (Any  other  numbers  that  may  be  assigned 
this  report) 


DISTRIBUTION    STATEMENT 


This  document  has  been  approved  for  public  release  and  sale;  its 
distribution  is  unlimited. 


SUPPLEMENTARY    NOTES 


12.    SPONSORING   MILITARY    ACTIVITY 

Naval  Postgraduate  School 
Monterey,  California  93940 


ABSTRAC  T 


The  background  of  microstructure  and  its  effects  on  sound  propagation  is 
reviewed.   The  expendable  bathythermograph  (XBT)  system  and  current  digitizing 
procedures  are  evaluated  to  determine  their  adequacy  in  investigating  vertical 
thermal  microstructure.   Several  methods  of  microthermal  analysis  are  proposed, 
Temperature  gradient  versus  depth  plots  are  used  to  analyze  data  from  several 
ocean  regions  -  Atlantic  Coastal,  Andaman  Sea,  Eastern  North  Pacific,  and  the 
northern  boundary  of  the  Gulf  Stream.   Analysis  of  these  data  suggest  a  water 
mass  dependence  of  the  thermal  microstructure.   The  significance  of  the  small- 
scale  structure  to  Naval  operations  is  discussed. 


DFORM 
1    NOV   68 

I    0101 -807-681 1 


1473 


(PAGE    1) 


Security  Classification 


A- 3 1408 


155 


Security  Classification 


key    wo  RDS 


Microstructure 
Microthermal  structure 

Expendable  bathythermograph 

XBT 

Digitizing 

Temperature 

Temperature  gradient 

Water  mass 

Underwater   sound 


DD  ,'M°o?..1473   back 


S/N    0101 -807-6H? 1 


156 


Security  Classification 


30  APR7I 

9  JUN72 

2  JUN  73 

7  JUL  73 

26  JUL  78 

18951 

S  1  0  1  3  7 

2  2  3  8  3 

2  12  9  0 

7  42 

Thesis 

120182 

M2765         Mag ruder 

c.1 

Some 

characteristics 

of 

temperature  micro- 

structure    in   the 

ocean . 

30APR71 

9  JUN72 

2  JUN  73 

7  JUI    7.-1 

1               26  JUL  78 

18951 

SI  0  1  3  7 

2  2  3  8  3, 
2  12  9  0 
2  4  7  42 

Thesis 

M2765 

c.l 


IEZ 


Mag ruder 

Somp  characteristics 
of    temperature  micro- 
structure    in   the 
ocean . 


thesM2765 

Some  characteristics  of  temperature  micr 

||||||||||||I|IMNII||1I!||I|!1II||PI|I 


3  2768  002  04185  7 

DUDLEY  KNOX  LIBRARY 


